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ABSTRACT 

We present the first detailed phase-resolved spectral analysis of a joint Chandra High Energy Trans- 
mission Grating Spectrometer and Rossi X-ray Timing Explorer observation of the p variability class 
in the microquasar GRS 1915+105. The p cycle displays a high- amplitude, double-peaked flare that 
recurs roughly every 50 s, and is sometimes referred to as the "heartbeat" oscillation. The spectral 
and timing properties of the oscillation are consistent with the radiation pressure instability and the 
evolution of a local Eddington limit in the inner disk. We exploit strong variations in the X-ray con- 
tinuum, iron emission lines, and the accretion disk wind to probe the accretion geometry over nearly 
six orders of magnitude in distance from the black hole. At small scales (1-10 R g ), we detect a burst of 
bremsstrahlung emission that appears to occur when a portion of the inner accretion disk evaporates 
due to radiation pressure. Jet activity, as inferred from the appearance of a short X-ray hard state, 
seems to be limited to times near minimum luminosity, with a duty cycle of ~ 10%. On larger scales 
(10 5 -10 6 i? g ) we use detailed photoionization arguments to track the relationship between the fast 
X-ray variability and the accretion disk wind. For the first time, we are able to show that changes 
in the broadband X-ray spectrum produce changes in the structure and density of the accretion disk 
wind on timescales as short as 5 seconds. These results clearly establish a causal link between the 
X-ray oscillations and the disk wind and therefore support the existence of a disk-jet- wind connec- 
tion. Furthermore, our analysis shows that the mass loss rate in the wind may be sufficient to cause 
long-term oscillations in the accretion rate, leading to state transitions in GRS 1915+105. 
Subject headings: accretion, accretion disks — black hole physics — instabilities — binaries: close — 
stars: winds, outflows — X-rays: individual (GRS 1915+105) 



1. INTRODUCTION 

Accreting stellar-mass black holes are known to exhibit 
different accretion 'states,' which are usually defined 
in terms of their X-ray spectr al shape and variability 
(jRemillard fe McClintockl 120061 and references therein). 
Physically, these states are intimately related to the fun- 
damental parameters of the accretion flow, e.g. the ac- 
cretion rate, the accretion disk geometry, the radiative 
efficiency of the disk, and the role of outflows in the form 
of winds and jets. The fact that X-ray binaries (XRBs) 
undergo state transitions highlights the dynamic nature 
of accretion onto black holes and allows for the possibil- 
ity of determining not only the fundamental parameters 
of accretion but also the physics that controls them. 

Of all the known Galactic black holes, GRS 1915+105 
is undoubtedly the most prolific source of state tran- 
sitions . Discovered as a transient by GRANAT in 
1992 ijCastro-Tirado et all ), it has remained in out- 
burst for the last 18 years and is typically one of the 
very brightest sources in the X-ray sky. It is also 
one of the most variable: its X-ray lightcurve con- 
sists of at least 14 different patterns of variability, 
most of which are h igh amplitude and highly-structured 
(iBellom et all 120001 he reafter BOO; iKlein-Wolt et al.l 
2002; Hamiik ainen et al.|[2005T) . It is believed that many 
of these phenomenologically-described variability classes, 
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which are labeled with Greek letters (|B00f ). are limit 
cycles of accretion and ejection in an unstable disk 
(IBellom et ail Il997at iMirabel et all 119981: iTagger et all 
2004: lFender fc Bellonill2004ft . 

The diverse timing and spectral properties of these 
variability classes make the behavior of GRS 1915+105 
particularly difficult to track in physical detail. Even in 
the x state, one of the simplest variability classes (with 
a relatively low X-ray flux, no structured variability, and 
a hard spectrum), GRS 1915+105 is never as faint or as 
hard as the ca nonical 'hard' state of black hole binaries 
(Belloni 2010). On the other hand, the canonical black 
hole 'hard' state is associated with jet production, and 
GRS 1915+105 has been shown to produce a jet during 
essentiall y every spectrally hard interval longer than 100 
seconds ( Klein- Wolt et al. 2002). This fact suggests that 
despite some differences in the spectral shape, the physics 
of jet form ation is probably the s ame for all Galactic 
black holes (|Fender fe Belloni|[200l . 

It is well-established that the properties of jets are 
highly correlated with the accretion state. But there 
is now a growing body of evidence suggesting that hot 
accretion disk winds are equally influenced by the accre- 
tion state, in part but n ot co mpletely due to chan ges in 
ionizing flux Lee et al.l 120021: iMiller et al.l l2006bl pOOl 
Neilsen fc Leell2009L hereafter Paper I, lUeda et al.ll2010l : 
Blum et al.ll2010D . In lPaperl . we demonstrated that the 
strength of the accretion disk wind in GRS 1915+105 is 
anticorrelated with the fractional hard X-ray flux, which 
is therefore a useful diagnostic of both the accretion state 
and outfl ow physics (for a pr eliminary discussion of these 
lines, see IMiller et al.1 12008). The anticorrelation holds 
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over many variability classes and there are indications 
that it may hold for o ther black holes a s well (XTE 
J16 50-500 fc GX 339- 4, IMiller et al.ll200l GRO J1655- 
40. IMiller et al J 12003: H1743-322. Miller et all 1200613 : 
IBhim et alJl2010h . 

If outflows from stellar-mass black holes depend on the 
accretion state, then any rapid variability should have 
observable consequences for those outflows. This is par- 
ticularly true for GRS 1915+105, where the X-ray spec- 
trum and the accretion flow can change drastically in sec- 
onds. For example, it has been shown conclusively that 
30-minute radio oscillations are 'baby jets' produced by 
ejecti o n events in cycles like the (3 state dFender et alj 
19971 iPoolev fc Fender] 119971 lEikenberrv et all 119981: 
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tion disk winds have been directly observed t o vary on 
simil ar (2-10 ks) timescales (GRS 1915+1 5: ILee et al " 



20021: lUeda. Yamaoka. fc Remillardl 120091 : lUeda et ah 
20Tq Cir X-l: iSchulz fc Brandt! I2002D or longer (1A 



0535+262: iRevnolds fc Miller! 120101 1 Flux-dependent 
stud ies have implied faster variations in both emission 
(1 s: IMiller fc Homanll2005l) and absorption lines (300 s; 
IMiller et al.ll2006bl ). However, the physical consequences 
of rapid variability on these winds have yet to be tracked 
in detail. A potential link between X-ray variability 
and disk winds would be especially interesting given our 
recent de monstrati on of a wind-jet interaction in GRS 
1915+105 (|Paper J ). 

In order to explore the relationship between jets, 
winds, and fast variations in the accretion disk, we have 
undertaken a detailed investigation of the p variability 
class in GRS 1915+105. Known affectionately as the 
'heartbeat' state because of the resemblance of its X-ray 
lightcurve to an electrocardiogram, the p state (see Fig. 
[T]) is a ~ 50-s oscillation consisting of a slow rise followed 
by a series of short bright bursts with amplitudes of 
order 5 x 10 38 ergs s _1 and strong changes in the 
X-ray spectral h ardness dTaam. Chen, fc Swank! 119971 
herea fter TCS97; iVilhu fc NevalainerJ Il998t iPaul et all 
119981) . Theoretical models suggest that this state is 
a manifestation of the Lightman-Eardley instability, a 
limit cycle in the radiation-pressure dominated inner 
accret ion di sk dLightman fc Eardlevl 119741: iBelloni et all 



1997al Uaniuk. Czernv. fc Siemiginowska 20001: 



Nava kshin. Rappaport. fc M elia 2000!: [janiuk fc Czernvl 
2005). 

Previously, we analyzed the average spectrum of a joint 
RXT E /Chand ra observation of GRS 1915+105 in the p 
state (Paper I). We found the average RXTE contin- 
uum during this observation to be relatively soft, with 
~ 79% of the 3-18 keV X-ray luminosity (L x ~ 4.9 x 10 38 
ergs s _1 ) emitted below 8.6 keV. In the time-averaged 
high-resolution X-ray spectrum from the Chandra High- 
Energy Transmission Grating Spectrometer (HETGS; 
ICanizares etHI l2005h . we detected an FeXXVl Lya ab- 
sorption line from the accretion disk wind with an equiv- 
alent width of -7.2 ± 1.7 eV and a blueshift of ~ 1400 
km s _1 . 

We follow up in this paper by tracking the p-phase- 
resolve timing and spectral variability of this same joint 
RXTE /Chandra HETGS observation of GRS 1915+105. 
For the very first time, we detect significant variations in 
absorption lines in phase-binned spectra, allowing us to 
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Fig. 1. — The 1-second PCA heartbeat lightcurve (top panel), 
HR1 (middle panel), and HR2 (bottom panel). HR1= B/A and 
HR2= C/B, where A, B, and C are the PCA count rates in the 
2.0-4.5 keV, 4.5-9.0 keV, and 9.0-30 keV bands. The heartbeat is 
composed of a slow rise or a shoulder followed by a double-peaked 
pulse. In HR1 the behavior is similar, although the shoulder follows 
the maximum rather than preceding it. 

assess real physical changes on timescales of seconds, well 
below the dynamical time in the wind. We explain the 
origin and evolution of the accretion disk wind via analy- 
sis of the X-ray spectral variability. Our results indicate 
that each bright burst has a significant impact on the 
accretion dynamics from the innermost to the outermost 
regions of the accretion disk. 

In Section [2] we describe our observations and data re- 
duction. In Section [21 we define the phase of the p cycle, 
explore variations in the recurrence time, and analyze 
phase-resolved power spectra. We perform joint spectral 
analysis with RXTE and the Chandra HETGS in Section 
[H We discuss our results on outflow formation and disk 
instabilities in Section (SJ summarize our understanding 
of the oscillation in Section [51 and conclude in Section [71 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra Data 

GRS 1915+105 was observed with the Chan- 
dra HETGS on 2001 May 23 (08:25:38 UT), for 30.16 ks. 
In order to mitigate pileup, the data were taken in Con- 
tinuous Clocking Mode; events were recorded in Graded 
format to reduce the risk of telemetry saturation. 

We reduced and barycenter-corrected the Chan- 
dra data using standard tools from the CIAO analysis 
suite, version 4.0. After reprocessing and filtering, we ex- 
tracted High-Energy Grating (HEG) spectra and created 
grating responses. We used the order-sorting routine to 
remove the ACIS S4 readout streak, since the destreak 
tool can introduce spectral artifacts for bright contin- 
uum sources like GRS 1915+105. We extracted 1-second 
lightcurves with dmextract. 

Due to incomplete calibration of Charge- Transfer 
Inefficiency (CTI) in CC/Graded mode, there is some 
wavelength-dependent disagreement in the continuum 
flux between spectral orders of the HEG (and the 
MEG, which we do not consider here because of its 
lower spectral resolution). For this reason, it is not 
currently possible to fit a physical continuum model to 
the HETGS data. Instead, we fit the individual spectra 
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Phase 

Fig. 2.— The phase-folded PCA heartbeat lightcurve in the (a) A = 2-4.5 keV, (b) B = 4.5-9.0 keV, and (c) C = 9.0-30 keV bands, 
along with (d) HR1 and (e) HR2, which are defined in Section [3] As discussed in the text, the mean period is about 50 seconds, so each 
major tick represents about 5 seconds. We have labeled the slow rise (~ 25 s), the soft and hard pulses, and the hard X-ray tail (~ 5 — 10 
s). Two cycles are shown for clarity. The gray points in the background are the individual data points, showing remarkable consistency 
despite noticeable variations in the cycle period. 



with polynomials to model the local continuum and 
use Gaussians for line features found in the combined 
residuals. 



2.2. RXTE Data 

During the Chandra observation, RXTE made a 
pointed observation of GRS 1915+105, beginning on 
2001 May 23 at 11:08:20.8 UT and lasting 21.9 ks 
(elapsed) with 13.7 ks exposure time. We select all avail- 
able data subject to the following constraints: (1) the 
Earth- limb elevation angle is above 3°; (2) the spacecraft 
is outside the South Atlantic Anomaly; (3) the offset an- 
gle from GRS 1915+105 is less than 0.02°. In this paper, 
we analyze the data from the Proportional Counter Ar- 
ray (PCA), which covers the 2-60 keV band. 

For timing analysis, we make use of the data from the 
binned mode B_8ms_16A_0_35_H_4P, which covers the 
2.0-14.8 keV band at 7.8 ms time resolution, and the 
event mode E_16us_16B_36_ls, which covers the 14.8- 
60 keV band at 15.3 fis time resolution. We extract 
1-second barycentered, dead-time-corrected, and back- 
ground subtracted lightcurves from each of these modes 
for timing analysis, and use our own software to cre- 



ate power spectra, subtr acting the dead-time-corrected 
Poisson noise level after iMorgan. Remillard. fc Greinerl 
(1997, hereafter MRG97). For high-S/N hardness ratios, 
we use the energy bands 2-4.5 keV, 4.5-9.0 keV, and 
9.0-30 keV, which have roughly equal count rates. 

For spectral analysis, the strong fast variability of the 
p state precludes use of the Standard-2 129-channel spec- 
tra, since the 16-second time resolution of the Standard-2 
data is too coarse for the 50-second cycles. To compen- 
sate, we use the binned mode and event mode data to 
create 32-channel spectra at 1-s time resolution (see, e.g. 
ITCS97L iBelloni et ai]ll997a[ ). We treat intervals with 
different combinations of PCUs separately. Our observa- 
tion features the PCUs operating in combinations {0, 2, 
3, 4} and {1, 2, 3, 4}, with roughly equal exposure times. 

3. TIMING ANALYSIS 

We show a representative portion of the PCA 1-second 
lightcurve in Figure [TJ along with PCA hardness ratios 
HR1 and HR2 (here defined as the ratio B/A and C/B, 
where A, B, and C are the count rates in the 2.0-4.5 
keV, 4.5-9 keV, and 9-30 keV bands, respectively). The 
lightcurve consists of a ~ 25 second slow rise, followed 
by a short burst, which is typically double-peaked (see 
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also Figure [5]). HR1 looks rather like the mirror image of 
the count rate (i.e. flipped horizontally about its maxi- 
mum) , except that HR1 is typically single-peaked. Other 
instances of the p state have shown different behavior, as 
the number and rela tive strengths of the peaks can vary 
(jMassaro et al.ll2010l Neilsen et al. 2011, in preparation). 

It is our main goal in this section to characterize 
the timing properties of the p state by tracking the ar- 
rival times of each individual burst. We define a phase 
ephemeris for the oscillation with <f> = at the time of 
maximum count rate (see the Appendix for the motiva- 
tion for this choice and the details of our method). This 
timing analysis is the foundation of all of the results pre- 
sented in this paper, as it allows us to study the evolution 
of the X-ray lightcurve, power spectrum, broadband X- 
ray spectrum, and the high-resolution X-ray spectrum as 
a function of phase (instead of time). 

3.1. Phase Ephemeris and Folded Lightcurve 

MRG97 tracked individual QPO waves in GRS 
1915+105 by fitting functional templates (e.g. sinusoids 
or Gaussians) to the PCA lightcurve. Given the high 
amplitude and unusual shape of the lightc urve desc ribed 
here, we opt for a modified version of the IMRG97I tech- 
nique (see the Appendix for more details). 

First, we take a single representative cycle from the 
data and cross-correlate it with the entire lightcurve. 
Maxima in the normalized cross-correlation values then 
indicate the times when = (peaks in the count rate). 
We then fold the data on this first set of (f> — times 
to obtain the average folded lightcurve. The process is 
iterated, with the folded lightcurve serving as a new tem- 
plate, to obtain the final set of <\> — times. This analysis 
results in 623 peak times from the Chandra lightcurve 
and 273 peak times from the PCA lightcurve, which 
should be accurate to 0.1 seconds or better for RXTE 
and 0.3 seconds or better for Chandra. The RXTE un- 
certainty is estimated from the variations produced by 
the use of different cross-correlation templates, while the 
Chandra uncertainty is derived from the offsets in <f> — 
times for cycles co-measured with RXTE (see the Ap- 
pendix for details). 

We also explored the 3-D pathway of the p cy- 
cle through the hardness-ha rdness-intensity diagram 
(HHID; see iSoleri et al.l 120081 for an example of this 
method) to consider using the mean pathway to con- 
struct a set of phase-resolved spectra that could be ana- 
lyzed and compared to the results derived from the cross- 
correlation ephemeris method. We find that the HHID 
method is less suitable for the p cycle because the sta- 
tistical noise for 1 s bins causes significant scatter in the 
HHID (particularly in the hard color dimension), and be- 
cause the rho cycle shows a tight loop in all three param- 
eters over the phase interval 0.04-0.12, causing some de- 
generacy for phase tracking based on the HHID. Thus we 
believe cross-correlations provide the optimal character- 
ization of the cyclic variability. With a phase ephemeris 
defined, we can create phase-folded lightcurves and hard- 
ness ratios (Figure [2|) . It is clear from this figure that 
phase-folding is a very effective way to characterize the 
oscillation, since the individual cycles have very similar 
shapes: the cycle-to-cycle variability is always less than 
20% for a given phase. 

To help diagnose whether there is systematic noise in- 
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Fig. 3. — The heartbeat period during our observations, mea- 
sured as the time difference between successive cycles. The peri- 
ods measured by Chandra (black) and RXTE (red) show very good 
agreement. Interestingly, the period exhibits considerable variabil- 
ity on short timescales and drifts slowly over long timcscalcs. 



troduced by using the broadband count rate to track 
the heartbeat cycles, compared to a parameter related 
to spectral shape, we conducted the following test. We 
used the PCA soft color, i.e. the ratio of source counts 
at 6-12 keV versus 2-6 keV, as an alternative quantity to 
compute cross correlations to define the times when each 
heartbeat cycle begins. The procedures were otherwise 
left intact. We note that the soft color usually shows 
one maximum per cycle, which is offset from the times 
of maximum count rate. The comparison of 239 cycle 
start times derived from these two tracking parameters 
(RXTE data; mean cycle period 50.488 s) yields an av- 
erage offset of 3.274 s, with a sample standard deviation 
of 0.283 s. It is the latter quantity that measures sys- 
tematic noise differences between the method that uses 
the count rate to track cycles, versus soft X-ray color. 
We conclude that neither our phase-resolved lightcurves 
and binned spectra nor our analysis conclusions would 
be significantly changed if the heartbeat phases were de- 
fined by soft color, since the phase deviations (0.0056) 
are much less than the spectral bin size (0.02). 

During our observation, each folded heartbeat is a slow 
rise (4> — 0.3 — 0.8 ~ 25 s), followed by a double-peaked 
pulse (</> = 0.8 — 0.1 « 15 s) and a hard X-ray tail lasting 
5-10 seconds (see <f> ~ 0.2 in panel c). Generally speak- 
ing, the slow rise is spectrally hard, the first pulse is soft, 
and the second pulse is moderately hard, but it is ob- 
vious even from Figure [5] that such statements must be 
made and interpreted only loosely, since HR1 and HR2 
are more or less a nticorrelated. 

Ot her authors (jJaniuk fc Czerny||2005t iMassaro et al.l 
120101 ) have already demonstrated that there is typically 
a delay between the soft and hard X-ray pulses, as well 
as a strong relationship between the times between these 
two bursts and the cycle period. Because HR1 peaks at 
the hard pulse, we see that the hard X-ray pulse lags the 
soft X-ray pulse by an average of 3.3 seconds in our ob- 
servation. This lag increases with the length of the cycle 
and appears to be a constant 6.5% of the cycle duration. 
This surprising result implies that the delay between soft 
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and hard X-ray pulses scales with the same clock as the 
entire oscillation. It also implies that using an X-ray 
color to define = simply introduces a phase shift 
into our results. Note that the variations in the lag are 
not large enough to interfere with our cross-correlation 
or phase-folding. We will return to the delay in Sections 
14.11 and 15.1.31 where we analyze the X-ray continuum. 

3.2. The Cycle Period 

We define the cycle period as a discrete quantity given 
by the time difference between two successive times of 
<j> = in the lightcurve. We show the heartbeat period 
for Chandra and RXTE in Figure [3] (black and red, re- 
spectively). Two points are immediately clear from this 
figure, both of which imply that the oscillation period 
in the p state is similar to that of low-frequency QPOs 
observed in GRS 1915+105: 

1. The heartbeat cycle period exhibits a slow secular 
drift from ~ 45 seconds to ~ 53 seconds over 30 ks. 

2. There is a substantial amount of scatter with an 
amplitude much higher than the secular drift. This 
scatter is also much larger than the uncertainty in 
the cycle period, which is of order 0.1 seconds for 
RXTE measurements (see Section l3~Tj) . 

The scatter in the period exemplifies the quasi-periodic 
nature of the heartbeat cycle: consecutive periods may 
vary by 5-10 seconds or more. The variability in Figure 
[3] is reminiscent of the behavior of QPOs in a variety 
of X-ray states of GRS 1915+105 as seen by IMRG971 
who tracked the arrival times of individual oscillations 
for several QPOs with frequencies from 0.067 Hz to 1.8 
Hz. 

MGR97 were able to demonstrate that the arrival 
times of these QPOs were well-described by uncorrelated 
Gaussian noise. For the p state, we find essentially no 
correlation between successive periods (the autocorrela- 
tion coefficient is r = — 0.018), and the cumulative distri- 
bution of periods is fit well with a two-Gaussian model. 
These Gaussians have mean periods p\ = 46.0 s and 
/it2 = 52.8 s, with standard deviations o~\ = 2.2 s and 

cr 2 = 2.1 S. 

The fractional scatter in the period measured here is 
compara ble to the p = 14.93 s and a = 0.72 s found by 
MRG97 for the 0.067 Hz QPO, and is actually less than 
the scatter measured for the other QPOs they tracked. 
On the other hand, the drift in the p period is much 
larger (roughly 40 phase units per 1000 cycles compared 
to 13 phase units per 1000 cycles for the 0.65 Hz QPO). 
Given its strong spectral variations, the p cycle is much 
more complex than a normal QPO, but the variability in 
its period is quantitatively similar. Given this similarity, 
it is interesting that a QPO at 6-15 Hz is present in the 
power spectrum during much of the p state, which is the 
subject of Section [331 

3.3. Power Spectra 

The phase-folded HR2 curve (Fig. 2) shows two max- 
ima that indicate enhanced hard X-ray components dur- 
ing the hard pulse ((f) = 0.05—0.1) and the hard X-ray tail 
(beginning near = 0.15). To help diagnose the spec- 
tral conditions during these intervals, we compute power 
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Fig. 4. — 2—10 Hz rms power as a function of cycle phase during 
the p state. We find rms > 0.1 (one of the criteria for the jet- 
producing X-ray hard state) in three phase windows: <f> = 0.02, rp = 
0.12—0.14, and in a broad interval near <j> = 0.3. While none of these 
windows corresponds to the hard pulse in the X-ray lightcurve, the 
broad interval does overlap with the hard X-ray tail. Thus the 
later stages of the hard X-ray tail may correspond to a short hard 
state. 

density spectra (PDS) for each second of source exposure 
with the RXTE PCA, using the count rate in the 2- 
37.9 keV band. In addition to subtracting the deadtime- 
corrected Poisson noise (see Section 2.2), we normalize 
the PDS to units of (rms / mean) 2 Hz -1 (|MRG97l ). 

To search for possible contributions from a jet during 
the heartbeat cycle, we hope to exploit the body of evi- 
dence that associates a steady radio jet with the proper- 
ties of the X-ray hard state. The hard state has the signa- 
tures of enhanced power (rms > 0.1) and a power contin- 
uum with "band-limited" shape, i.e., power density that 
appears flat at frequencies below a few Hz (TMunq et al.l 
l200lt lFenderll200a iRemillard k McClintockl 120061) . We 
note that measurements of rms in the literature often 
integrate the power over the range 0.1-10 Hz, but in this 
study a lower limit of 2 Hz is dictated by the Nyquist 
frequency for data samples taken at every second. 

We first examine the phase-folded values of rms power. 
This is derived by integrating each PDS over the range 
2-10 Hz, and then averaging the results in each of 50 
phase bins, using the heartbeat cycle phase ephemeris 
described in Section 3.1. The phase-folded rms curve is 
shown in Fig. [4] Here we see three maxima that exceed 
values of 0.1. These occur in narrow phase windows at 
0.02 and 0.12-0.14, and in a broad interval near <\> = 0.3. 
The first two maxima are likely associated with the in- 
trinsic variability of the oscillation, particularly the sharp 
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Fig. 5. — Power density spectra (PDS) corresponding to the three 
rms maxima in Figurel4l from top to bottom, <j> = 0.02, <j> = 0.12 — 
0.14, and ~ 0.26 — 0.36. During the first and second maxima, the 
PDS has a power-law shape. It is essentially featureless in the first 
maximum, but a QPO appears at 13 Hz in the second maximum. 
The third maximum shows band-limited noise that is consistent 
with the hard state, along with a QPO at ~ 8 Hz. 

decay of the soft and hard pulses, but we cannot rule out 
contributions from other sources of normal variability. 
However, none of the three intervals coincides with the 
hard pulse itself, which is instead aligned with a dip (to 
rms < 0.06) in Fig- [4] On the other hand, there is over- 
lap between the later portions of the hard tail and the 
broad interval of enhanced rms after <f> — 0.25. 

The PDS during the phases corresponding with the 
three rms maxima are shown in Fig. [3] Although the 
intrinsic variability of the p cycle may artificially en- 
hance the rms variability during the first two intervals 
and could bias the overall slope of the PDS, it is very un- 



likely to affect the qualitative shape of power continuum. 
At phase 0.02, the power density spectrum is featureless 
and the power density (P v ) declines with frequency (v) 
at a rate slightly steeper than P v oc v~ x . This PDS does 
not resemble the X-ray hard state. At the time of the 
second rms spike (phase 0.12-0.14), the power contin- 
uum again appears with a power-law shape, but there is 
additionally a QPO near 13 Hz. Finally, the PDS during 
the third phase interval with enhanced rms power does 
show a band-limited shape that resembles the hard state. 
Looking at the individual phase-binned PDSs, we find 
that similar shapes in the power continuum occur during 
the phase range 0.20-0.36, indicating the most likely in- 
terval for contributions from a jet that may temporarily 
form during each heartbeat cycle. 

The low-frequency QPO is detected during phase 
intervals 0.12-0.82 and 0.96-0.00, and it is therefore 
much more prevalent than either an elevated rms or a 
power continuum shape that resembles the hard state. 
The QPO frequency is 14 Hz at cj> = 0.12, after 
which it slowly decreases to 7.3 Hz near <f> = 0.52, 
and ends near 8 Hz at <j> = 0.82. As suggested 
by the bottom panel of Figure [5j a first harmonic is 
usually present. The presence of a QPO while the 
source evolves through a long, hard dip in the X-ray 
lightcurve is reminiscent of the (3 cycle in GRS 1915+105, 
which is a 30-min cycle tied to the formation of im- 
pulsive radio jets dMarkwardt. Swank, fc Taaml 119991: 
Mikles Eikenberrv. fc Rothsteinl 120061: iMirabel et"all 
1998t iFender fc Bellonil I2004D . Given that the X-ray 
spectra during the /3 cycle were successfully interpreted 
as the combination of thermal radiation from an accre- 
tion disk plus a hard X-ray power law attributed to in- 
verse Compton emission, the QPO behavior in the long 
hard dip of the p cycle may justify application of the 
inverse Compton model to this state. 

However, this leaves behind two puzzles. First, the 
hard pulse coincides with a phase interval that exhibits 
neither high rms values nor a low-frequency X-ray QPO, 
and so we find no temporal signatures that might pro- 
mote spectral interpretation of this hard X-ray compo- 
nent via Comptonization. Secondly, there is another 
brief appearance of a QPO at phases 0.96-0.00 (with fre- 
quency that shifts from 6 to 7 Hz). This QPO is weaker 
than the 7-14 Hz feature and lacks a first harmonic, so it 
is difficul t to link this short interval to th e longer QPO 
episode. iSoleri. Belloni. fc Casellal (|2008[ ) also detected 
transient QPOs in GRS 1915+105, associating them with 
state transitions. But because our transient QPO does 
not coincide with an obvious state transition, the relation 
between these features is unclear. 

Finally, in some phase intervals (0.02—0.06, 0.46—0.50, 
and 0.68 — 0.74) there is also weak evidence (2.8 — 3.1c) 
of a QPO (rms ~ 1%) near 60 Hz. This is of inter- 
est be cause GR S 1915+105 is known to exhibit a 67 Hz 
QPO (MRG97) that can be particularly strong in 7-type 
lightcurves. A hint of this feature is seen at 59 Hz in the 
top panel of Fig. [51 At other phases the feature looks 
more like an edge, as suggested by the bottom panel of 
Fig. [5] Examination of the PDS in individual phase bins 
suggests that this edge is not an artifact of a stronger 
QPO moving to lower frequency. These features are too 
weak to support further analysis, and so this topic is left 
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Fig. 6. — PCA spectra and residuals for two different phases of 
the heartbeat cycle. The fainter spectrum from the cycle minimum 
(phase = 0.3) is significantly fl atte r than the brighter spectrum 
at the peak, = 0. See Section 14.11 for details of our continuum 
models. 

to a future generation of instruments with larger collect- 
ing area. 

4. SPECTRAL VARIABILITY ANALYSIS 

In the following subsections, we explore in detail the 
spectral variability of the p state, using RXTE to study 
the broadband X-ray properties of the oscillation and 
the Chandra HETGS to probe the known accretion disk 
wind as a function of p cycle phase. In the previous 
section, we de monstrated fas t spectral evolution in the p 
state (see also TCS97; BOO), but it remains to be seen 
if and how accretion disk winds and jets participate in 
this variability. Here we will show that the changes in 
the X-ray continuum are related to significant variations 
in the accretion disk wind, both occurring on timescales 
< 5 seconds. This constitutes the very first probe of disk 
wind physics on s uch short timescales. All spectral fittin g 
is done in ISIS (IHouck fc Denicolal f2000l IHouck! l200l . 
We assume a distance and incli nation of D = 11.2 kpc 
and i = 66° (|Fender et al.l[l999h : we fix N H = 5 x 10 22 
cm~ 2 ()Lee et al.l 120021 and references therein). 

4.1. RXTE PCA 

For our variability analysis, we apply our derived phase 
ephemeris from Section 13.11 to extract phase-resolved 
spectra: for each 1-second PCA spectrum, we com- 
pute the p cycle phase and average the results using 
50 phase bins. At each phase, we fit the spectrum 
from 3.3 to 45 keV. X-ray spectra of black hole bina- 
ries have been modeled as the combination of various 
soft and hard components, which may include but are 
not limited to: thermal radiation from the accretion 
disk, power law emission, an explicit treatment of Comp- 
tonization, or bremsstrahlung emission. For the highly 
variable states of GRS 1915+105, a model consisting of 
a multi-temperature accretion disk plus some type of 
power law has been especially effective at tracking the rel- 
ative ch anges in the accretion disk and Comptonizat ion 
(see e.g. IBelloni et ai]|1997at IMigliari fc Belloni|[200l . 

We focus here on the variations of the accretion disk 
and hard X-ray components inferred from several differ- 
ent continuum models. We show two example phase- 
resolved PCA spectra (from the peak and minimum of 



the oscillation) in Figure|51 As expected from our energy- 
resolved lightcurves, the spectrum at the peak of the 
count rate (0 = 0) is significantly brighter and softer 
than the flat spectrum at the cycle minimum (<f> = 0.3). 

4.1.1. Model 1: simpl 

Our initial model for the X-ray continuum con- 
sists of five component s: cold absorption (tbabs; 

iWilms. Allen, fc McCravl I2000D. a hot disk component 
(ezdiskbb: iZimmerman et al.ll2005l ). an emission line at 
6.4 keV (egauss), a high-energy cutoff (highecut, whose 
formula is given by exp(E/ E^u ) ) , plus a scattering com- 
ponent (simpl: TSteiner et al.ll2009l ). simpl is a convolu- 
tion model that takes any seed spectrum and scatters 
a fraction /sc of the photons into a power law. The 
high-energy cutoff is required to account for curvature 
in the hard X-ray spectrum, and we choose simpl be- 
cause it conserves photons and avoids the divergence of 
the powerlaw model from realistic expectations of Comp- 
tonization at low energies. We allow the different PCU 
combinations to have different parameters, although we 
tie the high-energy cutoff parameters together for bet- 
ter constraints. In the end, our best fit parameters for 
the different PCU sets typically differ by < la. Model 
1 provides an excellent statistical description of the X- 
ray continuum at nearly all phases of the p cycle (overall 
X 2 /^ = 2301.9/2150 = 1.07), but the parameters during 
the hard pulse (4> ~ 0.06) are somewhat puzzling. See 
below and Section 14.1.31 for more details and alternative 
explanations. 

We show the results of our broadband spectral fits in 
Figure [7] In these fits, the observed maximum temper- 
ature in the disk T ^ s hovers around 1.1-1.2 keV for the 
majority of the cycle, spiking sharply to ~ 2.2 keV during 
the X-ray burst. At the same time, the inner radius of 
the disk is relatively large (70-110 km) and grows dur- 
ing the slow rise, dropping swiftly to ~ 40 km at the 
end of the soft pulse. The inner radi us R\„ is related to 
the e zdiskbb normalization -/V disk as ([Zimmerman et aLl 
2005): 



Adisk — -Jl 



D 



(1) 



Because simpl counts photons when calculating the disk 
normalization and the power law (/sc)i Adisk in Model 
1 uniquely tracks the disk parameters prior to Comp- 
tonization. For the high luminosities here, we use a color 
correction factor / = 1.9 (J. Steiner, private communi- 
cation). The resu lts are qualitatively s imilar to the fast 
variations seen bv IBelloni et al.l (jl997blla? ) in the n and A 
states, although their disk radii are a factor of ~ 2 lower 
than ours. This difference appears to be mainly due to 
our use of ezdiskbb rather than diskbb. Thus it does 
not indicate a physical difference between the p state and 
the k and A states. 

During the slow rise the photon index V is stable 
around 2.1, which is relatively steep for a hard state 
but not unusual for GRS 1915+105. But near = 0, 
the r displays a double-peaked pha se dependence, risin g 
quickly to its upper limit (r = 4; iSteiner et al.l 12009). 
dropping sharply back to ~ 2.5 during the hard pulse, 
then peaking again near 4. After the pulses, T decays ex- 
ponentially back to 2 with an e-folding interval A<fi ~ 0.1. 
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Fig. 7. — Model 1 fit parameters and la error bars as a function of phase for the PCA spectra (averaged over different PCU combinations). 
Two cycles are shown for clarity; since the average period is 50.33 seconds, each phase bin corresponds to roughly 1 second of real-time 
variability. The te mperature and radiu s of the accretion disk vary significantly throughout the cycle, reminiscent of changes in the ft/A 
states analyzed by Bclloni ct al. (ll99Tbl faT> . Because simpl counts disk photons to create the power law ( /sc)i Model 1 uniquely tracks the 
disk parameters prior to Comptonization, i.e. scaling up the disk normalization by a factor 1 — /so See 14.1 II for details of the model fits. 



Meanwhile, the scattering fraction /sc decreases slowly 
during the slow rise and reaches a minimum ~ 0.13 near 
<j> = 0.95. Then, during the hard pulse at <j) ~ 0.06, 
/sc spikes sharply to 1; there is a second, broader peak 
with /sc ~ 0.6 during the hard X-ray tail (Fig. [7]). We 
interpret this behavior as an indication that during the 
hard pulse, the X-ray spectrum is completely dominated 
by scattering and there is essentially no direct disk com- 
ponent. Still, the disk dominates the light during the 
soft pulse, contributing as much as 70% of the observed 
flux just before = 0, and roughly 25% during the hard 
X-ray tail. 

The high-energy cutoff exhibits similarly strong varia- 
tions around the phase of the hard pulse. During most 
of the cycle, -Efoid is steady around 30 keV, although the 
spectrum of the hard X-ray tail is effectively consistent 
with no cutoff, or with -Efdd < 1 MeV. But during the 
hard pulse, £fdd drops sharply to 6 keV. Thus Model 1 
leads to the interpretation of the hard pulse as a period of 
very strong (i.e. Compton thick) scattering by relatively 
cool electrons. 

The flux (i*6.4 kov) m the Gaussian emission line (see 
Section I5.2[) declines very slowly for most of the cycle, 



but peaks during the soft pulse. With the exception of 
this short pulse, the line flux is very strongly correlated 
with the line width <J6.4 keV, which varies between about 
0.5 keV and 1.5 keV. The correlation coefficient for these 
two parameters is r — 0.88. At the RXTE spectral res- 
olution, it is difficult to reliably decouple the variations 
in the line width and the line flux, so the true variabil- 
ity of the feature is unclear. The line is too broad to be 
observable in the Chandra HETGS spectrum. However, 
if we include only the phases of the cycle away from the 
pulses, 4> = 0.15 — 0.8, we find a moderate correlation 
between the line flux and the simpl scattering fraction 
(r ~ 0.63). We will return to this in Section I5~2l 

4.1.2. Model 2: nthcomp 

In the interest of quantifying the robustness of 
our results to the choice of hard X-ray component, 
we also explore the evolution of the disk with sev- 
eral prescriptio ns for Comptonization, including compTT 
and nthcomp (jZdziarski. Johnson, fc Magdziarzl 119961 : 
iZvcki. Done, fc Smith| [l999i Here we discuss the results 
from our modeling with nthcomp. To use this model, we 
set the seed photon distribution to be a disk blackbody 
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Fig. 8. — Model 2 fit parameters and la error bars as a function of phase for the PCA spectra (averaged over different PCU combinations). 
Two cycles are shown for clarity; since the average period is 50.33 seconds, each phase bin corresponds to roughly 1 second of real-time 
variability. The disk temperature and radius are very similar to those of Model 1. The nthcomp seed photon temperature is tied t o the 
accretion disk temperature, so that both parameters are well-constrained even when the disk normalization goes to zero. See Section l4.1.2l 
for the details of Model 2. 



and tie the seed temperature to the disk temperature. 
This allows the disk temperature to be well-constrained 
even when the disk normalization goes to zero (see be- 
low) . The resulting variations are qualitatively very sim- 
ilar to those measured with simpl, and are shown in Fig- 
ure^ The overall goodness of fit is xt — 1-09. The power 
law index T and the electron temperature kT e are strik- 
ingly similar to the corresponding parameters in Model 
1, and although the normalization of nthcomp is defined 
differently than /gc, these two parameters also behave 
alike. See Table [1] for a comparison of the fit parameters 
at 4> = 0.06. 

Still, there are several noticeable distinctions between 
Model 1 and Model 2, which are primarily related to 
the fact that with simpl, the photons in the hard X-ray 
component come directly from the disk component, while 
nthcomp treats them as distinct. For example, we find in 
Model 2 that our disk radii and temperatures are some- 
what lower, peaking at approximately 100 km and 1.8 
keV instead of 110 km and 2.2 keV. Furthermore, during 
the hard pulse, the spectrum is so completely dominated 
by the nthcomp component that the disk normalization 
(and radius) go to zero. This phenomenon is certainly 



unphysical, so we ultimately prefer simpl (and its use of 
fsc t° track the pre-Comptonizcd disk) to nthcomp as a 
description of the spectrum. 

Nevertheless, the strong changes in the accretion disk 
and Comptonization parameters around the hard pulse 
appear to lead to the same conclusion as Model 1. In 
both models, the hard pulse is characterized by strong 
scattering, a relatively steep photon index (r = 2.5 in 
Model 1, 3.2 in Model 2), and a population of scatter- 
ing electrons with significantly-reduced energy (-Efoia = 6 
keV in Model 1, kT e = 4-5 keV in Model 2). These num- 
bers are not characteristic of the canonical X-ray hard 
state. Therefore, while it is clear that Comptonization 
models do provide a good description of the X-ray spec- 
trum, our findings here confirm our conclusion from Sec- 
tion 13.31 that it is difficult to interpret the hard pulse as 
a brief X-ray hard state. 

4.1.3. Model 3: Bremsstrahlung 

We are now left with a small dilemma. The Comp- 
tonized disk models provide a very good statistical de- 
scription of the phase-resolved X-ray spectra of the heart- 
beat state. However, several points indicate the possi- 
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TABLE 1 

X-ray Continuum Properties at <f> = 0.06 (Hard Pulse) 



Parameter 


Model 1 


Model 2 


Model 3 


Rin (km) 


41 ± 2 


< 12.4 


22.1 ± 0.4 


J- obs (Ke V ) 


7 194. n no 
L. 1Z IE U.UO 


1 7 A -4- n 09 
1 . t 41 31 u.uz 


z. ry m u.uo 


r 


2.86+°'" 

ww — 0.07 


3.29+°'°3 

u ^-0.29 


> 3.2 


£ 

JSC 


> 0.99 




n nnv+0 27 
0.007_ 004 


-^ntiiconip 




6.23+°'!" 

—0.76 




■^bremss 






9.56 ±0.06 


S fo i d (keV) 


8-4-0.9 






kT c 




4 - ' u -0.07 




^brcmss (keV) 






f- no + 0.05 


^6.4 keV 


027+ 008 

u - uzl -0.009 


012+ 006 


0.051 ±0.008 


06.4 keV (keV) 


< 1.0 


< 1.2 


0.8 ±0.2 


x 2 /^ 


21.9/43 


24.5/43 


32.0/40 



Note. — The hard pulse is the point of the p cycle that has 
the most unusual spectral and timing properties, and coincides 
with the most dramatic changes in the continuum parameters. 
It therefore provides an ideal contrast between our models. The 
values reported here are averages of the two PCU combinations; 
the individual errors (90% confidence ranges for a single parame- 
ter) have been added in quadrature with the standard deviation 
of the values for the two PCU sets. % 2 is calculated for 56 data 
points. 

ble presence of a second hard component in the X-ray 
spectrum. The most powerful argument in favor of such 
a component comes from the X-ray spectra themselves. 
For more than 70% of the cycle, both the photon index 
and the electron temperature are completely constant, 
and the scattering fraction varies smoothly. But these 
parameters change dramatically during the hard pulse 
and then return to "normal." This suggests that there 
may be two geometrically-distinct sources of hard X-rays 
(one roughly constant and one variable). Additional evi- 
dence for a second component comes from the power den- 
sity spectra fSection l3.3j) : during the hard X-ray tail and 
the beginning of the slow rise, the rms level, the shape of 
the PDS, and the evolution of the QPO frequency are all 
consistent with the properties of the hard-state-like dips 
in GRS 1915±105. These dips are known to be domi- 
nated by Compton scattering. The hard pulse (Fig.[2p), 
on the other hand, exhibits a dip in rms and there is no 
low-frequency QPO, which is typically associated with 
Compton-dominated intervals at high luminosity. This 
further indicates that the mechanism producing the hard 
X-rays during the hard pulse may be distinct from that 
in the hard X-ray tail and the slow rise. Finally, our 
timing analysis of 242 RXTE observations of the p state 
(Neilsen e t al. 2011, in preparatio n) a nd the theoretical 
model s of iNavaksh in et al.1 (|2000( ) and Uaniuk fc Czernyl 
(|2005t) support the presence of a second component dur- 
ing the hard pulse that may be related to plasma ejec- 
tions from the inner disk. 

To test this idea, we replace the high-energy cutoff in 
Model 1 with a thermal bremsstrahlung component to 
represent the emission from this ejected plasma. The 
overall quality of the fit is comparable to the other mod- 
els (xl = 1-12), but because simpl alone provides a 
decent fit outside the X-ray pulses, the hard compo- 
nent is occasionally over-determined. Nevertheless, this 
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Fig. 9. — Luminosities of spectral components as functions of 
cycle phase. In the top panel, for Model 1, the bolometric disk 
luminosity (before scattering) is shown as a dot-dashed black line; 
we plot the 3-45 keV total, observed disk, and scattered (power- 
law) luminosities as red, green, and blue solid lines, respectively. 
The middle panel is the same but for Model 2, with the luminosity 
of the nthcomp component in blue. The bottom panel (Model 3) 
additionally shows the bremsstrahlung component in purple. 

model presents a physically-interesting alternative to the 
Compton-dominated models, so we discuss it here briefly. 

The accretion disk temperature and inner radius be- 
have similarly to the parameters shown in Figure [71 
although we find again that the disk normalization is 
smaller due to the presence of the bremsstrahlung com- 
ponent. But in light of the sudden changes in the hard 
X-ray components during the hard pulse in Models 1 
and 2, we are primarily interested in the phase depen- 
dence of the bremsstrahlung normalization and temper- 
ature in Model 3. Here, we find a sharp spike in the 
bremsstrahlung normalization from K = 1 — 2 at = 
to a value of K = 9.56 ± 0.06 at cj> = 0.06. The normal- 
ization is given by 



K 



3.02 x 10" 15 
4nD 2 



(2) 



where D is the distance to GRS 1915±015, n e and rii 
are the electron and ion number densities, and V is the 
emitting volume. The bremsstrahlung temperature dur- 
ing this spike is fcTbrcmss ^5 — 7 keV, which is consistent 
with the evolution of the electron temperature in Models 
1 and 2. Although the exact physical interpretation of 
the hard X-ray component varies from model to model, 
all three models require the sudden appearance of a new 
population of electrons during the hard pulse, which coin- 
cides with dramatic changes in the accretion disk. These 
electrons are cooler than typical coronal electrons, but 
warmer than the disk. 

For comparison, we show the luminosities of the X-ray 
continuum components in Figure from which several 
points are immediately clear. First, the hard X-ray lu- 
minosity is relatively steady in all three models, with 
the exception of the hard pulse, which is scattering dom- 
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Fig. 10. — Accretion disk bolometric luminosity versus tempera- 
ture for our continuum models (Model 1 in black, Model 2 in red, 
Model 3 in green) . We use a star to mark the location of <j> = for 
both curves; motion around the loop is clockwise. We also over- 
plot dashed and dot-dashed blue lines of L oc T 4 and T 4 / 3 , which 
are expected if the disk radius or the accretion rate is constant, 
respectively. Note that we have cut off Model 2 (nthcomp) where 
the disk luminosity goes to zero. 

inated in Models 1 and 2 but 56% bremsstrahlung in 
Model 3. Second, the hard X-ray tail (</> = 0.15 — 0.3) is 
also dominated by scattering (here the disk component 
constitutes < 30% of the light in all models). Finally, we 
can see that in all models, just before = the bolo- 
metric disk luminosity (see also 14. 1.41 for details) reaches 
80-90% of the Eddington luminosity for GRS 1915+105 
(a 14 M Q black hole; L Edd = 1.8 x 10 39 ergs s" 1 ). Model 
1 actually tops out at 1.03LEdd- 

4.1.4. Disk Variations and the Accretion Rate 

The most obvious feature of our fits to the heartbeat 
state X-ray continuum is the strong variability in the 
temperature and inner radius of the accretion disk. Be- 
fore turning to a theoretical interpretation of the p cycle 
(Section 5), we assess the implications of our continuum 
fits for the mass accretion rate onto the black hole. 

Given the color correction factor / and the measured 
inner radius and temperature of the disk, we can calcu- 
late the disk mass accre tion rate using equation (4) of 
Zimmer man et ail (|2005l) : 

/ • \ 1/4 

obs ~ 2.05 ^ 8TrR? n a J ' W 

where G is the gravitational constant, M is the black 
hole mass, and a is the Stefan-Boltzmann constant. We 
evaluate the accretion rate as a comparative scale, not- 
ing that the absolute value is quite uncertain, and that 
we assume no changes in spectral hardening or the ra- 
diative efficiency of the disk. All three models give very 
similar mass accretion rates and bolometric luminosities. 
We find that M peaks just before the spike in the disk 
temperature at <fi = 0, at a value M < 1.5 x 10 19 g s _1 . 
After the peak, the mass accretion rate drops sharply to 
M ~ (1 - 3) x 10 18 g s" 1 . For all models the evolution of 
M is similar to that of the bolometric disk luminosities 
(Figure [9]). 
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Fig. 11. — HETGS model/data spectra for three different phases 
of the heartbeat cycle with the rest energies of Fe xxv (Hca: Is 2 - 
ls2p, 6.7 keV) and Fexxvi (Lya: Is - 2p, 6.97 keV) marked in red 
and black, respectively. At (f> = 0.2, near the cycle minimum, the 
iron absorption lines are weak. At <j> = 0.54, we detect noticeable 
features from Fexxv and Fexxvi, both blueshifted by ~ 500 km 
s _1 . The feature at 6.2 keV may be Mnxxiv Heo or some other 
Doppler-shifted absorber. At cj> = 0.92, the Fexxv line is gone, but 
there is some evidence for faster (~ 5000 km s — 1 ) and broader iron 
absorption in the Fexxvi line profile. 

The evolution of the accretion disk component can also 
be seen in Figure 1 101 where we plot its bolometric lu- 
minosity versus its temperature. During the p cycle, 
GRS 1915+105 traces a clockwise loop in the L — T 
plot. In Models 2 and 3, the loop is distorted be- 
cause of the competition between the disk, nthcomp, and 
bremsstrahlung components, but the behavior is other- 
wise si milar. Although we calcu late each point on the 
plot ala I Zimmerman et al.l (|2005l ): 

L = 73.9a f^V i?? n = 73.9 aT 4 bs ^^, (4) 
V / / cos * 

there is essentially no interval longer than a few seconds 
over which L oc T 4 . After <j> = 0, which is marked with a 
star, the disk moves rapidly from its moderate tempera- 
ture, high luminosity state through a high temperature, 
moderate luminosity state {<fi ~ 0.06) to a low tempera- 
ture, low luminosity state (<fi ~ 0.2). After the minimum, 
the luminosity rises at roughly constant temperature (in 
our nomenclature, this interval is the slow rise; see Fig. 
(2b). For reference, we overplot lines of L oc T 4 , which 
results from a constant inner disk radius, and L oc T 4 / 3 , 
which is expected for a constant accretion rate. These 
lines broadly match the evolution of L going into and 
out of the slow rise. We will return to this plot, its in- 
terpretation, and the question of the mass accretion rate 
in Section [5j 

4.2. Chandra HETGS 

While the RXTE PC A provides high S/N broadband 
X-ray spectra for measurements of the continuum, the 
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Fig. 12. — (a) Measured flux (in units of 10 -3 photons s" 1 cm~ 2 ) 
in the Fe XXV (red) and Fe XXVI (black) absorption lines as a func- 
tion of p cycle phase. The H-like iron line clearly tracks the X-ray 
continuum, although the peak line flux follows the peak X-ray lu- 
minosity with a delay Acj> = 0.92 fa 46 seconds. In contrast, the 
He-like line only appears briefly, growing in strength with Fexxvi 
until <f> ~ 0.54, where it starts to fade, (b) The full band RXTE 
PCA phase-folded lightcurve, for ease of tracking changes in the 
absorption lines observed with Chandra (top panel). 

Chandra HETGS provides an excellent high-spectral- 
resolution characterization of the observed narrow lines 
in the soft X-ray band that is ideal for studying the ac- 
cretion disk wind and its dynamical evolution. 

One of the most significant challenges for fast phase- 
resolved grating spectroscopy is choosing time and phase 
intervals to maximize the detectability of interesting fea- 
tures and their possible variations. For ease of compari- 
son with the PCA spectra, we restrict our attention here 
to the ~ 20 ks of the Chandra observation that bracket 
the RXTE data. This choice maximizes the S/N in the 
FeXXV line, which increases over the course of our ob- 
servation, relative to Fexxvi (see Section l5.3.4l) . For our 
phase-dependent spectral analysis, we extract 50 HEG 
spectra evenly spaced in phase. In order to achieve suf- 
ficient S/N, we set the phase width of each spectrum to 
A(f> = 0.2. This is essentially a sliding box window, so 
there is so me overlap between cons ecutive spectra, but as 
shown by iSchulz fe Brandt] (|2002f) , this method is quite 
suitable for line variability studies. As discussed in Sec- 
tion 12.11 because of calibration uncertainties we model 
the Chandra X-ray continuum with a polynomial fit (in- 
stead of the physical models of Section l4~T]) . Polynomials 
accurately characterize the local continuum, so they ef- 
fectively isolate narrow absorption lines. See Figure [TT] 
for three example residual spectra. 

A complete search for variable features at all veloci- 
ties and ionization levels with any arbitrary phase de- 
pendence is beyond the statistical scope of this paper. 
Since we are mainly intereste d in the known accretio n 
disk wind in GRS 1915+105 (ILee et alJl200l IPaper ID . 
we focus on the observed narrow iron absorption lines 
with velocities ranging from ~ 500 to 2000 km s _1 (see, 
e.g. the middle panel of Figure [TT]). We detect both 
Fexxvi (Lya: Is - 2p, 6.97 keV) and Fexxv (Hea: Is 2 - 
ls2p, 6.7 keV). We use Gaussian fits to measure the pa- 



rameters of these iron absorption lines at each phase (see 
Table[2]). To decouple our results from possible variations 
in the wind dynamics, which could manifest as changing 
line widths and thus affect the apparent line flux, we fix 
the intrinsic line width at 200 km s _1 . This is a typical 
orbital speed in the ou ter disk, and matc hes the turbu- 
lent line width found by Ue da et al.l ()2009f ) in their study 
of the 4> state. When we detect both Fe xxv and Fe xxvi, 
we can fit for a 90% confidence upper limit on the line 
width of 800 km s -1 . 

Finally, we note the possible presence of additional 
lines that vary with phase, including a line at 6.2 keV 
that may be Mnxxiv Hea and an unidentified feature 
at 4.29 keV. Without secure identifications it is unclear 
if they are part of the iron absorber or if they represent 
a different dynamical component, such as an accelera- 
tion zone of the wind or a signature of infall. Both have 
maximum fluxes around 2.0 ±0.5 photons s _1 cm -2 ; the 
corresponding equivalent widths are 2.5 eV and 6.4 eV, 
respectively. Since we cannot clearly identify their ori- 
gin we do not consider them further here, beyond noting 
that their phase dependence is somewhat similar to that 
of the iron absorber. 

It is clear from both Table [2] and Figure [TTJ that the 
ionization parameter £ of the wind changes substantially 
during the cycle. Tabled] shows that Fexxv and Fexxvi 
peak with similar equivalent widths at very different 
phases ((f) = 0.54 vs. <f> = 0.92), while a comparison of the 
middle and bottom panels of Figure [TTJ shows that the 
Fexxvi absorption line is strong even when the Fexxv 
line has disappeared. In the rest of this section, we fo- 
cus on these two later phases as representative of the 
variations in disk wind absorption in this state. 

This ionization evolution is even more obvious in Fig- 
ure [121 The Fe XXV absorption line appears to track the 
Fexxvi line in the interval = 0.2 - 0.54, but it fades 
rapidly outside this interval and is not significantly de- 
tected for the rest of the cycle. The Fexxvi absorption 
line, however, continues to grow until <j> = 0.92. The peak 
flux in this line is more than twice that measured from 
the average spectrum; Monte Carlo simulations indicate 
that this increase is significant at 96% confidence. 

It is also evident in Figure [T2] that the Fe XXVI absorp- 
tion line evolves similarly to the PCA count rate. Under 
the assumption that the line variability is a response to 
the X-ray continuum impinging on the wind, it appears 
that the peak line flux follows the X-ray flux with a delay 
A(j> w 0.92 = 46 seconds. This is actually rather long: 
the semima.jor axis of GRS 1915+10 5 is a ~ 250 lt-s 
(|Greiner. Cubv. fe~ VIcCaughrcan 200lj), and previous es- 
timates ha ve placed the wind near r ~ 10 l t-s from the 
black hole (|Lee et al.ll200llUeda et al.ll2009UPaper J ). It 
may then be reasonable to suspect that A<f> is mainly due 
to continued increases in the ionization of the absorbing 
gas. We will show in Section [531 that changes in the gas 
density are also important. 

Additionally, we note that our measurements of the 
absorption line centroids provide marginal evidence for 
variations in the blueshift of the wind during the cycle. 
At <p = 0.54, we detect Fexxvi at v ~ -AZV^q km s -1 ; 
by cj) = 0.92, the velocity has increased to v ~ —1600^50° 
km s _1 . This apparent acceleration is significant at 96% 
confidence, and it occurs in combination with the ap- 
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TABLE 2 

X-ray Absorption Lines in the p State of GRS 1915+105 



Phase 


Line 


E 
(keV) 


E hs 

(keV) 


(km s" 1 ) 


Flux 


W 
(eV) 


0.20 
0.54 
0.92 


Fexxvi Lya 
Fe xxv Hea 
Fe XXVI Lya 
Fe xxv Hect 
Fe xxvi Lya: 
Fe xxv Hea 


6.966 
6.700 
6.966 
6.700 
6.966 
6.700 


7 01+ - 02 
'• ul -o.oi 

6.711 a 

6 976+ 018 

' ±1 -o.oo8 

7 no4+ - 001 

6.711 a 


_ 17 00+soo 

-490 a 
-430+ 390 
-490+ 370 
-1610±f™ 

-490 a 


< 1.6 

< 0.5 
1.5 ±0.5 

i r+0.4 

3.9 ±0.7 

< 0.1 


> -8.6 

-0 1+ ' 1 
u - 1 -2.6 

— 1 1 o +3 ' 5 

— 10 1+ 2 - 2 
±u ' ± -3.0 

-14 1+ 2 ' 5 
> -0.5 



Note. — Errors quoted arc 68% confidence ranges for a single parameter. Phase: line 
parameters are reported for three different phases of the cycle; Eo : rest energy; E ^ s : 
measured energy; AiJ s hif t : measured Dopplcr velocity; Flux: measured absorbed line flux in 
units of 10 — 3 photons s~ 1 cm~ 2 ; Wq : line equivalent width; max : phase of maximum line 
flux. Boldface indicates the phase of maximum flux for that ion. 

a Since the Fcxxv line is not detected at this phase, we fix its energy at the value measured 
at <j> — 0.54 to place limits on the line flux and equivalent width. 



pearance of a blue wing to the Fe xxvi profile (see Fig. 
[TTj) . We then take the velocity and line width as some 
evidence of variability in the dynamics of the disk wind. 
In Section [5.31 we discuss the photoionization evolution, 
dynamics, and the 46-second delay in the context of wind 
formation scenarios and the disk-wind-jet connection. 



5. DISCUSSION 

In this section, we explore the significance of our re- 
sults from Sections [3] and |U Specifically, we will address 
our findings of (1) strong, fast variations in the inner ra- 
dius of the accretion disk, (2) bremsstrahlung emission 
during the hard pulse, and (3) clear changes in the ab- 
sorption lines from the accretion disk wind on timescales 
of 5 seconds. 

While strong variations in the X-ray spectrum of GRS 
1915+105 are routinely observed, the magnitude and 
short timescale of the spectral variability in the heart- 
beat state bears repeating, especially given its cyclic na- 
ture. In the rising phase of the cycle, the bolometric disk 
luminosity changes from about (5 to 16) x 10 38 ergs s _1 
in roughly 35 seconds; during the decay phase, the X-ray 
luminosity decreases by as much as Lx < 5 x 10 38 ergs 
s~ 2 . In other words, we find that the black hole halts a 
burst reaching ~ 80 — 90% of its Eddington luminosity in 
just a few seconds. As discussed in Section 4, the lumi- 
nosity variations are associated with substantial changes 
in the inner radius and temperature of the accretion disk. 

The physical changes in this bizarre oscillation are rep- 
resented in the phase variation of the X-ray continuum, 
fluorescent emission lines, and the hot iron absorption 
lines of the accretion disk wind. In the following subsec- 
tions, we discuss our results on phase-resolved spectra 
and power spectra to understand the origin, dynamics, 
and accretion geometry of the p state and to explore in 
detail the apparent coupling between the accretion disk 
wind and the X-ray oscillation. We describe our model 
for the physics of the heartbeat state, starting with varia- 
tions in the accretion rate in the inner disk, and tracing 
their influence on the accretion dynamics all the way to 
the outer edge of the disk. 



5.1. X-ray Continuum 
5.1.1. M Variations and Z Source-like Behavior 

The most obvious properties of the heartbeat state 
(based on our X-ray continuum fits) are the observed 
strong variability in the luminosity, temperature, and in- 
ner radius of the accretion disk (Figs.[7 HTU|) . Specifically, 
our analysis shows that during the slow rise, the disk ra- 
dius gradually grows, then drops sharply and rebounds 
during the X-ray peaks. Meanwhile, the disk tempera- 
t ure appears constant and t hen spikes rapidly. 

ILin. Remillard. fe Homanl (2009, hereafter LRH09) 
discovered similar variations in the radius, temperature, 
and luminosity of the accretion disk in the accreting neu- 
tron star XTE J 170 1-462 as it evolved from a Z-source 
to an atoll source. Z sources are accreting neutron stars 
that display three branches in the shape of a "Z" in their 
color-color diagrams. LRH09 found that in the ensemble 
of observations in the vertex between the flaring branch 
and normal branch, the inner radius of the disk increased 
with luminosity above L ~ 0.2LEdd, while the disk tem- 
perature remained roughly constant. They interpreted 
this result as an indication of a local Eddington limit 
(IFukuell2004tlHeinzeller fe D uschl 20Q3) in the inner disk. 
Local Eddington effects arise because the radiation forces 
and gravitational forces depend differently on radius in 
thin disks, so that there exists a critical radius where 
radiation pressure can truncate the disk. Inside this crit- 
ical radius, a significant fraction of the disk is expelled by 
radiation pressure. The resulting spectrum is the same 
as that of a standard thin disk for energies E > fcsT cr i t . 
It is possible that the same process is at work during the 
slow rise of the p cycle in GRS 1915+105, i.e. when the 
radius increases with luminosity at roughly constant T 

(Fig. HOD. 

LRH09 also showed that on the XTE J1701-462 flaring 
branch, the luminosity is proportional to T 4 / 3 , which is 
expected if the inner radius of the disk changes while the 
accretion rate remains constant. They concluded that 
the flaring branch is an instability in which the disk tem- 
porarily slips inwards but is quickly driven back to the 
equilibrium position seen at the vertex position on that 
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same day. In their interpretation, this equilibrium posi- 
tion is set by the local Eddington limit. For comparison, 
we plot a line of L oc T 4 / 3 in Figure [TO] It can be seen 
that in GRS 1915+105, the disk very roughly approxi- 
mates this scaling as it both approaches and exits the 
slow rise of the p cycle. We also show L oc T 4 ; the two 
lines bracket the behavior of the accretion disk during 
the p cycle. These rough scalings seem to indicate that 
some changes in the disk may occur at ^constant accre- 
tion rate, especially after <j> = 0. 

Thus, mapping LRH09's prescription for neutron stars 
onto GRS 1915+105, we can understand the gradual in- 
crease in the disk radius during the slow rise of the heart- 
beat state as a local Eddington phenomenon (our mea- 
sured radii and tempera tures are consistent with those 
predicted byjFukuc 2004J). The disk slowly expands with 
an increasing accretion rate. After the accretion rate 
reaches a maximum the disk first falls in at ~ constant 
accretion rate and then cools rapidly while the accretion 
rate drops, setting the stage for the next cycle. It should 
be noted that the analogy between the XTE J1701-462 
lower vertex/flaring branch and the p cycle is not exact, 
since the variations in GRS 1915+105 are much faster, 
and XTE J 170 1-462 does not trace out a loop in the 
L — T diagram. Furthermore, only GRS 1915+105 shows 
the rapid catastrophic changes seen in the disk and the 
accretion rate around = 0. 

Nevertheless, it is very exciting to see such similar 
accretion disk behavior between a neutron star system 
and a black hole system; the commonalities suggest that 
the p cycle may share or mimic some properties of a 
Z source. This particular comparison seems reasonable, 
since Z sources are known for high accretion rates (see 
LRH09 and references therein), and GRS 1915+105 may 
have the highest accretion rate of any known black hole 
transie nt. These similarities also support t he conclu- 
sion of lDone. Wardzihski. fc Gierlihskil (pOOl that GRS 
1915+105 may be unique merely for its high accretion 
rate. 

5.1.2. On the Role of Radiation Pressure 

Our interpretation that the local Eddington limit in 
the inner accretion disk may produce the slow rise 
during the p cycle (see Section 15. 1 . 1|) poses an in- 
triguing question as to the role of radiation pres- 
sure in the heartbeat state. We now have claims 
of two very different radiation pressure phenomena in 
the p state: a local Eddington limit and (histori- 
cally) the Lightman-Eardley radiation pressure instabil- 
ity (RPI). The RPI is commonly invoked, with some 
success, to expl ain the strong limit cycle oscil lations in 
GRS 191 5+105 (lTCS97t lNavakshin et al.ll2000l hereafter 
NRM00: [Janiuk fc Czernvl l2005L hereafter JC05). These 
processes are physically independent but not necessarily 
conflicting. Here we consider the question: what is the 
true role of radiation pressure in the p state? Does it 
produce limit-cycle oscillations by making the disk (in- 
tact at constant Ri n ) thermally/ viscously unstable, or 
does it actually push out the inner radius of the disk 
over timescales of ~ 25 s w 4000 dynamical times? Or 
are some elements of both models require to explain the 
observations? To answer this question, we need to un- 
derstand how the observed changes in the radius, tem- 
perature, and luminosity of the accretion disk might be 



related to these physical mechanisms. 

For the local Eddington limit, R[ n increases with L at 
constant T. For half of each observed p cycle (the slow 
rise; Fig. EJd), the disk luminosity and radius appear to 
grow at ^constant T. Alone, however, a local Eddington 
limit cannot produce an oscillation or a loop in the L — T 
diagram. Some other mechanism is required to cyclically 
vary the accretion rate for a complete description of the 
heartbeat state. 

The question is more subtle as it relates to the RPI. 
Does the RPI act in concert with a local Eddington limit 
to produce the p state, or can it mimic an interval of 
increasing i?j n at constant T? Since time-dependent disk 
simulations of the RPI (e.g. NRM00 and JC05) have not 
produced synthetic spectra, the question is also difficult 
to answer. For preliminary considerations, we examine 
the dynamics of the accretion disk in their simulations. 

Physically, the RPI occurs when the inner disk be- 
comes unstable due to radiation pressure. This drives a 
limit cycle involving an oscillation of the local accretion 
rate in the inner disk. This can be interpreted as a 'den- 
sity wave', which originates around i?dw = 20 — 30i? g . 
Here R g is the gravitational radius (21 km for GRS 
1915+105). As the density wave moves through the disk, 
it partially evacuates a gap in the disk. Inside and out- 
side the gap, the disk surface density is enhanced (see 
Figure 4 in both NRM00 and JC05). The inner density 
enhancement peaks around R = 5R g ~ 105 km. The 
similarity of this estimate to our measured radii is strik- 
ing. During the decay of the cycle, this excess density 
moves inwards quickly (and our measured radii decrease 
sharply). These results provide a hint as to why the ap- 
parent disk radius might change rapidly in the pulses. 

It is unclear, however, if the RPI can mimic our obser- 
vation of increasing i?; n at constant T during the slow 
rise. Our local Eddington interpretation implies that 
the true edge of the disk moves outwards with L, while 
NRM00 and JC05 have R in fixed at 2 R g . On the other 
hand, the disk surface density and temperature are in- 
creasing and decreasing functions of radius, respectively, 
and both groups find a short hot state and an extended 
"cold" state. For these reasons, it seems potentially pos- 
sible for their time-dependent thin disk models to mimic 
an expanding-radius disk at constant T. 

Thus the role of radiation pressure in the heartbeat 
state requires some additional model analyses. A local 
Eddington limit in the disk does a good job explaining 
the slow rise, but cannot explain the pulses or the ex- 
istence of an oscillation. The RPI is a viable explana- 
tion for the oscillation, but it is unclear if it can either 
act alone or jointly with the local Eddington limit to re- 
produce our observations of the slow rise. Future time- 
dependent disk simulations will be able to resolve this 
issue by calculating model spectra that can be compared 
with the observations. For now, it seems that our ob- 
served disk radii and temperatures provide some support 
for both the RPI and a local Eddington limit at work in 
GRS 1915+105. 

5.1.3. Ejections, Jets, and the Fate of the Inner Disk 

Another exciting question posed by the observed rapid 
oscillations in the radius, temperature, and mass accre- 
tion rate in the disk concerns the ultimate fate of the 
excess material in the inner disk. Does it simply fall into 
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the black hole, or is it ejected from the inner disk and 
therefore potentially observable? We have not directly 
observed any ejection events, but the mysterious hard 
pulse in the X-ray lightcurve still lacks a satisfying in- 
terpretation. It is even peculiar theoretically, since RPI 
simulations produce only single-peaked cycles. Is this 
short burst of hard X-rays related to the ejection of ma- 
terial from the inner disk, or does it have an explanation 
in terms of standard black hole states? 

The association of the hard pulse with a typical soft- 
to-hard state transition appears to be ruled out, since 
the spectral and timing properties of the hard pulse are 
not remotely consistent with typical hard states for ei- 
ther GRS 1915+105 or black hole binaries in general. 
That is, although we are able to successfully model its 
X-ray spectrum with disk/Comptonization models, the 
resulting parameters are inconsistent with the canoni- 
cal hard state, and imply the sudden appearance of a 
low-temperature, high-optical-depth cloud of electrons 
(t ~ 5 with nthcomp). 

Before going ahead, we may ask if the appearance of 
these cooler electrons can be explained simply by an 
increase in their cooling rate. If synchrotron cooling 
is important, additional cooling could be accomplished 
by an increase in the electron density or magnetic field 
strength. If inverse Compton cooling dominates, then 
an increase in cooling could result from the observed in- 
crease in the X-ray flux or from an increase in the electron 
density. Detailed radiative transfer calculations would be 
necessary to determine the energetics precisely, but it is 
unclear in the Compton cooling case why the cooling is 
is abrupt, while the flux changes smoothly. The simplest 
explanation in both cases is a sudden increase in the elec- 
tron density. This justifies our use of a bremsstrahlung 
component to model the hard pulse, and provides us 
with an observational basis for discussing plasma ejec- 
tions from the inner disk. 

Additional support for plasma ejections comes from 
the theoretical work of NRM00 and JC05, who found 
that the X-ray states of GRS 1915+105 are reproduced 
with higher fidelity when ejection processes are included. 
NRM00 allowed a fraction of the total accretion power 
to be channeled into a jet with Lorentz factor 7^3. 
JC05, modeling the p state specifically, allowed the disk 
to evaporate into the corona at the peak of the cycle. 

It should be n oted that these models (see also 
Uaniuk et al.l 12000) produce plasma ejections when the 
X-ray luminosity is gr eatest, and are therefor e at odds 
with the conclusions of Kle in- Wolt et al.l (|2002T ) . who ar- 
gued that jet production is essentially a continuous pro- 
cess that occurs du ring hard faint intervals (State C in 
the classification of IBOOI ). However, the p state is only 
known to exhibit low-level radio emission, and it is un- 
clear if this is a 'baby' jet with a short duty cycle or some 
other stable but faint radio act ivity. Thus there may be 
no observational conflict with Klein- W olt et al.l (2002) 
for the p state. Future high time-resolution radio/X-ray 
observations (e.g. with the EVLA and Chandra) should 
probe the precise nature of the radio emission in such 
states. 

In any case, both groups predict plasma ejections near 
= 0, and our models implies the sudden appearance 
of some cold plasma around the same phase. With 
the reasonable assumption that these two points are 



related, we proceed to estimate how much material is 
ejected from the disk. For this, we can use equation ([2]), 
which relates the bremsstrahlung normalization to the 
volume and electron density in the emitter. Because the 
bremsstrahlung pulse begins just after the drop in R ln 
at <fi ~ 0, we assume that i?t>rcmss ~ max(i?i n ) ~ 97 km, 
so that the emitting volume is V meiK ~ 4 x 10 21 cm 3 . 
In other words, we assume that some bremsstrahlung- 
emitting material is ejected from the inner disk, filling a 
sphere of radius 97 km. Assuming the particle density 
is constant over the sphere and that the material is full y 
ionized (with ISM abundances from I Wilms et al.l [2000). 
we find an implied ion density ri\ ~ 10 20 cm -3 and ion 
mass m\ ~ 10 18 g. For a Shakura-Sunyacv disk with 
M ~ 10 19 g s _1 , the disk mass interior to 97 km is of 
order 10 21 g, > 1000 times larger than the mass of this 
bremsstrahlung cloud. Since JC05 argued that as much 
as 10% of the disk accretion rate may evaporate into the 
corona, our requirement that 0.1% of the disk mass con- 
tribute to the bremsstrahlung pulse seems reasonable. 

Also of interest is the possible production of a colli- 
mated jet at some point in the p cycle. In Section 13.31 
we considered the extent to which the X-ray minimum 
of the p state, or even the slow rise, could be considered 
a hard state, with the purpose of identifying those phase 
intervals where we might expect to see some signatures of 
jet production. Our examination of the phase-dependent 
PDS and X-ray continuum indicates that the end of the 
hard X-ray tail, <j) = 0.26 — 0.36, has timing and spec- 
tral properties that are consistent with an X-ray hard 
state. Thus if there is any jet activity, we might expect 
it to have a duty cycle of order 10%. In this window 
around the count rate minimum, a hard state (and jet) 
may temporarily coexist with the disk. This conclusion is 
particularly interesting in light of the possibility of strong 
synchrotron cooling during the hard pulse, but a radio 
detection of a synchrotron spectrum would be required 
to verify and connect these two phenomena. 

Previously, we reported the Ryle telescope radio 
flux near this observa t ion to be Sisghz 5 mJy 
(|Poolev fc Fenderl Il997t iNeilsen fc Lei [2009) . which is 
comparable to the sensitivity of the Ryle telescope. Thus 
we cannot draw any robust conclusions about a jet. We 
note, however, that this is averaged over nearly a full cy- 
cle; if the actual radio flux is close to this upper limit, 
then the "true" jet strength could be an order of mag- 
nitude higher. Thi s would be typical for a \ state jet 
in GRS 1915+105 dFender et al.(fl999t iKkdn-Wolt et al.1 
2002). Since a short-lived jet is plausible, future sensi- 
tive radio observations of the p state may place strong 
constraints on jet evolution around black holes. 

To summarize, it appears that our Comptonized disk 
plus bremsstrahlung model leads to a physically-sensible 
observational interpretation of the p state in the context 
of plasma ejections. We find that the double-peaked X- 
ray lightcurve is the result of two processes: (1) large, 
cyclic variations in the mass accretion rate in the inner 
disk, leading to an increase in the disk temperature (the 
soft pulse), and (2) the ejection/evaporation of a small 
fraction of the inner disk, which produces a brief burst 
of bremsstrahlung as it propagates into the corona (the 
hard pulse). The hard X-ray emission could also sim- 
ply be scattering in a Compton-thick cloud, but the ap- 
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Phase 

Fig. 13. — Flux (top) and equivalent width (bottom) of the broad 
iron emission line as a function of phase. The line flux is roughly 
constant during most of the cycle and spikes during the two pulses, 
while the equivalent width oscillates slowly. 

pearance of this cloud is still best explained by material 
evaporating from the disk. This "steam" from the disk 
does not qualify as a jet, but it is nevertheless an exciting 
signature of plasma ejections during the heartbeat cycle. 

5.2. Iron Emission Line 

Here we briefly consider the broad iron line seen in 
our RXTE spectra (Figure [T3J) . In our models, the 
line flux slowly decreases for most of the cycle and the 
equivalent width oscillates slowly, spiking during the soft 
pulse. Models 1 and 2 show a dip in the line strength 
during the hard pulse, but the reality of t his feature is 
unclear. In the standard interpre t ation dFabian et al 
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20031 and references therein), the broad line is a reflec- 
tion feature from the inner accretion disk. As shown in 
Section 14. 1[ the innermost regions of the disk fluctuate 
cyclically during the p state, so the fact that the line flux 
is almost constant indicates that some geometrical effects 
from the oscillation may be influencing the illumination 
of the disk. For example, the increase of the disk scale 
height with luminosity may shield outer regions of the 
inner disk from some of the X-ray luminosity, allowing 
the reflected line flux to stay the same. 

On the other hand, if the line originates at the inner 
edge of the disk, then the constant line flux might be 
caused by a combination of increased illumination and 
decreased reflecting area. We note that the line equiva- 
lent width looks rather like the simpl scattering fraction, 
so it seems likely that the changes in the hard X-ray com- 
ponent are literally reflected in the iron line. But without 
detailed theoretical models of the disk spectrum and the 
iron line profile, a robust determination of its location is 
not possible. Thus it would be difficult to use this broad 
line to place a constrain on the black hole spin. In other 
states like the x state, where the X-ray variability and 
its effects on the disk structure are not so important, 
estimates of the spin paramete r a* may be both more 
feasible and mor e meaningful (|Martocchia et all 120021 : 
IBlum et al.l[2009h . 



5.3. Accretion Disk Wind 

One of the main goals of this work, as a follow-up to 
our discovery of a wind-jet interaction in Paper 1, is to 
understand the complex relationship between the fast X- 
ray variability and the physics of the accretion disk wind, 
i.e. its origin, dynamics, ionization, structure, and short- 
timescale evolution. This includes determining whether 
or not the structure of the wind is constant in time. One 
could easily imagine a scenario in which the wind is un- 
affected by the limit cycles in the inner disk, e.g. if it 
is launched from the outer disk by simple X-ray heating 
from the phase-ave raged X-ray luminosity o r by MHD 
processes (see, e.g. IMiller et all I2006al H008). Or per- 
haps there is a simple relationship between the wind and 
the variability, where the wind forms independently but, 
as it rises off the disk, is photoionized by the strong quasi- 
periodic X-ray variability. 

In this section, we exploit the simultaneous strong vari- 
ability in the accretion disk wind (Figure I12[) and the 
X-ray continuum to argue that neither of these scenar- 
ios provides a satisfactory description of the accretion 
disk wind physics in GRS 1915+105. We use photoion- 
ization arguments to probe the accretion dynamics and 
plasma conditions in the outer disk, and conclude that 
there must be a powerful coupling between the formation 
of accretion disk wind and the X-ray luminosity. We ar- 
gue that photoionization alone is insufficient to produce 
the observed variability in the wind. 

5.3.1. Plasma Conditions and Ionization Balance 



In Section 14.21 we showed that the absorbed fluxes in 
both Fexxv and Fexxvi vary with </>. Given the vari- 
ability in the X-ray absorber, we treat the disk wind 
as though it is not (in general) in photoionization equi- 
librium. However, it is reasonable to assume that the 
observed maxima in the absorbed fluxes occur when the 
instantaneous photoionization and recombination rates 
are momentarily equal. Here, under this assumption, we 
estimate the gas density in the wind. 

Neglecting collisional ionizatio n, this tempor ary ion- 
ization balance can be expressed (|Liedahllll999D 
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(5) 



where R is the distance to the continuum source, n c , ni, 
and are the electron density and number density of 
iron in charge states i and i + 1, a^+i is the recombina- 
tion rate into charge state i (out of charge state i + 1), 
and $i is the ionizing spectrum integrated over the pho- 
toionization cross-section o~i for charge state i : 
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Here L £ is the monochromatic luminosity and \i is 
the ionization threshold. Given a set of ion fractions 
(or Hi and n i+ i), an ionizing continuum, and appro- 
priate recombinati on rates, the ionization parameter 
(jTarter et aU[l969[) 
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is determined by equations (|5|) and 
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In general, the ion fractions depend on the ionizing con- 
tinuum, the ionization para meter, and the overall density 
(|Kallman fe Ba utista 200l|) . However, if the absorption 
lines are unsaturated, we have an additional constraint 
from the linear part of the curve of growth, which re- 
lates the equivalent width of a spectral line to the ionic 
column density: 

W\ ire 2 

= ^NiXfr. (8) 
A m e c z 

Here W\ is the line equivalent width in mA, A is the 
wavelength, m e is the electron mass, c is the speed of 
light, Ni is the column density of charge state i, and 
is the oscillator strength of the relevant transition. 

At 4> = 0.54, the Fexxvi line is almost certainly un- 
saturated because it continues to increase in both ab- 
sorbed flux and equivalent width for another 20 seconds 
(for saturated lines, the flux and equivalent width grow 
very slowly until the optical depth is > 1000). For this 
reason, it seems likely that the Fe xxv line is also unsat- 
urated. Using our measured equiv alent widths and os- 
cillat o r strengths from the XSTAR (|Kallman fc Bautista! 
120011: iBautista fe Kallmanl 1200 ID database, we evalu- 
ate equation [8] for Fexxv and Fexxvi. We find that 
nxxvi/«xxv « Nxxvi/Nxxv = 2.5 at <p = 0.54. 

To recap, under the reasonable assumptions that the 
wind is optically thin and in temporary ionization equi- 
librium, we can estimate the density by specifying the 
location of the wind (n e oc R~ 2 ). Based on the dynam- 
ical estimates of IPaper 1 . we set i? w ind ~ 2 5 x 10 11 cm. 
We us e photoionization cross-sections from IVerner et al.l 
(1996) and recombination rates from xstar (T. Kall- 
man, private communication). For temperatures 10 K 
< T < 10 7 K, we find that the implied electron density in 
units of cm -3 is 10.8 < log(n c ) < 12.4. For the measured 
luminosities we can expect the ionization parameter to 
be between 10 3 and 10 5 . This is reasonable for a plasma 
dominated by hydrogen- and helium-like iron. 

5.3.2. Ionization Parameter and Plasma Dynamics 

At the beginning of Section lSTBl we posed the question: 
is it possible that the structure of the wind is constant 
in time, so that the variability in the spectral lines is 
purely due to the changing ionizing luminosity? Here we 
show that the answer is a resounding 'no'. This implies 
that the wind must be re-formed or re-launched each and 
every cycle. Even without the photoionization consider- 
ations that follow, this answer could be anticipated from 
Figure where it can be seen that between = 0.2 
and (j) = 0-54, the flux in the Fexxvi and Fexxv lines 
are roughly equal, while the X-ray luminosity increases 
by ~ 40% over the same phase interval. It is quite un- 
likely that these observed line fluxes could be produced 
purely by photoionization, because the ion fractions for 
H- and He-lik e species are only compara ble over narrow 
intervals in £ ()Kallman fe Bautistall200lT ). 

That said, we performed a brief photoionization anal- 
ysis with XSTAR, using the analytic model warmabs. 
Since we are essentially modeling only H- and He-like 
iron, it should be understood that these models are de- 
signed not to give an exhaustive account of the photoion- 
ization state of the gas but to parametrize the relation- 
ship between the X-ray continuum and the wind. We use 



our models of the ionizing spectrum (Section l4.ip at <f> = 
0.54 and <j) — 0.92 to generate atomic level populations. 
We then fit our 6-7.5 keV HETGS spectra for the equiva- 
lent hydrogen column density, ionization parameter, and 
velocity shift of the wind at those phases. We fix the 
turbulent line width at 200 km s _1 (Section 14.21 see also 
lUeda et al.ll2009D . We set the density at n = 10 12 cm" 3 . 
Our best fits indicate that at (j) = 0.54, when Fexxv 
is maximized, the equivalent hydrogen column density, 
ionization parameter, and velocity of the accretion disk 
wind are A H = 2.1±g;£ x 10 22 cm~ 2 , log£ = 3.87^;™ 
ergs cm s _1 , and v = —550^390 km s -1 . The errors 
are la. At <fi = 0.92, the maximum of Fexxvi, we find 
Ah = 10±l x 10 22 cm" 2 , log£ = 4.8 ±0.1 ergs cm s -1 , 

and v = -1520^240 km s_1 for tnis disk wind - More 
detailed analysis would be necessary to account for ad- 
ditional dynamical components at this phase (Fig. Ill|). 

We have already pointed out the changing velocity of 
the absorber, but we notice here a significant increase 
in the column density and ionization parameter of the 
wind (factors of roughly 5 and 7, respectively). During 
this same phase interval (<fi = 0.54 — 0.92) the ionization 
parameter increases far more than the X-ray luminosity. 
The luminosity increases by a factor of 2, and the ioniz- 
ing flux for H- and He-like iron only increases by a factor 
~ 1.4 — 1.5. This is an indication that the luminosity vari- 
ations are insufficient to produce the observed absorption 
line variability. The apparent changes in the fitted col- 
umn density are further evidence that the structure of 
the wind cannot be constant in time. This, too, could 
be anticipated from Figure II H where both the wind ion- 
ization state and the strength of the Fe xxvi line change 
significantly with (f>. 

Admittedly, the exact value of the column density 
could be underestimated if resonance scattering or ther- 
mal emission lines fill in some of the a bsorption features 
durin g the heartbeat cyc le (see, e.g. iWoidowski et al.l 
120031: iKallman et alj[2009h . Additionally, the column of 
Fe XXVI may be enhanced by recombination in the totally 
ionized wind. However, even if these effects are impor- 
tant, significant ionization evolution in the disk wind is 
still implied by the fact that Fexxv disappears while 
Fexxvi remains. The appearance of a blue wing in the 
Fexxvi line profile near <f> — 0.92 (Figure [TTj) is a final 
indication that the structure of the outflow varies sub- 
stantially during the cycle, probably due to the injection 
of additional material. 

5.3.3. Wind Formation Scenarios 

The observed strong variability in the accretion disk 
wind leads to two intriguing conclusions. First, the den- 
sity and structure of the wind must change on very short 
timescales. This must be the case because the wind be- 
comes progressively more ionized as the cycle proceeds, 
but the X-ray luminosity does not increase enough to 
produce the resulting ionization. Similarly, the distance 
between the wind and the X-ray source cannot change 
much in the requisite ~ 20 seconds because the outflow 
is too slow. Thus equation ([7]) dictates that the den- 
sity or structure of the wind must change. Second, these 
changes in the wind must be quasi-periodic (as observed). 
Fast quasi-periodic variability in the wind is rather un- 
usual, so we address it first. There are three basic meth- 
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ods by which the wind might show coherent variability 
during the cycle: 

1. Its structure is constant in time, but its ionization 
parameter changes in accordance with Lx- 

2. The wind is periodically launched from within or 
close to the oscillating region of the inner disk. 

3. The wind is periodically launched from the outer 
disk (where previous studies have suggested it orig- 
inates) by radiation from the inner disk. 

Based on simple ionization arguments, we have already 
ruled out Case 1. Cases 2 and 3 are more complex sce- 
narios, and we weigh their pros and cons in what follows. 
We are unable to rigorously rule out these scenarios, but 
the outer-disk origin seems preferable (see below). 

For Case 2, let us suppose the wind originates at 50 R g . 
Since the phase variations in this region are strong, an 
oscillating wind seems natural. However, given our mea- 
sured ionization parameters, the density in the wind 
would have to be n c ~ 3 x (10 17 - 10 19 ) cm" 3 . Neither 
thermal nor radiation-driven winds can produce such 
dense, highly ioni z ed ou tflows at such small radii (see 
iProga fc Kallmanl (|2002l ) and references th e rein), which 
would leave MHD processes (|Progal 120001: iMiller et all 
\m)6A\2(M\ . But even if the magnetic field does vary 
coherently during the p cycle, it is difficult to understand 
how a wind at 50i? g could have a velocity as low as 500 
km s" 1 and a line width < 800 km s _1 , given that the 
orbital velocity at 50 R g is approximately 40,000 km s" 1 
for a 14 Mm black hole with a spin parameter a* = 0.98 
(jMcClintock et afl l2006) . It is also unclear how this wind 
could persist for 30 ks given the short infall time at 50 R g . 

For Case 3, let us suppose that the wind originates 
at 10 lt-s from the black hole. At this radius, dynami- 
cal and viscous timescales are much too long to produce 
coherent phase variability. However, the strong contin- 
uum variability provides a straightforward origin for the 
wind: as the bright X-ray burst propagates outwards, it 
should irradiate the disk in ever-larger annuli, launching 
a thermally driven wind from each radius. This imprints 
a strong phase dependence on the wind. In other words, 
the outer disk and inner disk are linked by radiation. 

The most significant challenge for this scenario involves 
getting the wind into the line of sight quickly, since the 
inclination of the disk (i.e. the jet axis) is 66°. There 
are several options to resolve this dilemma, all of which 
amount to reducing the angle between the plane of the 
outer disk and the line of sight: (a) Irradiated disks are 
known to flare vertically at large radii, with scale heights 
H/R < 0.1; (b) radiation- or tidally-driven wa rps may 
be present in the outer disk of GRS 1915+105 (jPringld 
119961: iWhitehurst fc Kiiig 119911 ): (c) the inne r disk and 
outer disk may be misaligned, as suggested bv lLee et al.l 
(2002) to connect the overabundance of iron along the 
line of sight to the geometry of the sup ernova that formed 
the black hole (e.g iBrown et al.l l2000). With a combina- 
tion of these possibilities, it may be possible to produce 
the observed 30-40 second rise time for the wind. How- 
ever, it should be noted that this scenario, while favored, 
cannot be fully tested at this time because our HETGS 
phase resolution is limited. We find that shifting our 



= times by some constant offset (on top of the ran- 
dom scatter reported in Section I3.2[) merely introduces a 
phase offset into Figure [12] Longer observations of this 
wind may be able to measure the travel time directly. 

Case 3 also has the attractive potential to explain 
qualitatively a number of the observed properties of 
the disk wind, especially its increasing column density, 
ionization parameter, and velocity. The column den- 
sity increase is simple: more material is launched into 
the light of sight as time passes. In a thermally-driven 
wind, the velocity and ionization parameter are propor- 
tional to R 1 , and the density decre ases like n oc i?~ 2 ~ 7 
(|Begelman. McKee. &: Shields! 1 1983} ) . 7 is a power- law 
index between and 0.5 under equilibrium conditions. 
If such scaling laws are applicable in our observation, 
then the increase in v and £ may be partly related to 
thermal driving of the wind. A robust estimate of 7 
would require us to determine the precise launch time of 
the wind, or the point at which the optical depth in the 
newly-launched wind dominates that of the very highly- 
ionized wind from the previous cycle. Such a probe is 
beyond the capabilities of our present data. Simulations 
of the oscillation-driven wind (Proga et al., in prepara- 
tion) may shed additional light on these issues. 

5.3.4. The Role of the Wind in GRS 1915+105 

Finally, we attempt to estimate the mass-loss rate 
in the wind in our observation. This is particularly 
important because rece nt hydrodynamic si mulations of 
thermally-driven winds ([Luketic et al.|[20Toh have shown 
that such winds can have mass loss rates well above the 
accretion rate and thus may exert a strong influence on 
the accretion dynamics of the system. Conditions here 
are likely to be far from steady state, so we use an el- 
ementary consideration based on the maximum column 
density in the wind. In the cylindrical approximation, 
which is appropriate if the wind is launched ^vertically 
off the disk, M w ind should scale approximately as 

M w i nd =74rnA^H/At launch (9) 
= &nR w - mA v z mN- a , (10) 

where Aii aunc h is the time over which the wind is 
launched, v z is the wind speed perpendicular to the plane 
of the disk, A — 2 x 2irR w i n dV z At\ aviIiC } l is the outer sur- 
face area of the cylinder filled by the wind, and m is the 
average ion mass per hydrogen atom (~ 2.4 x 10 ~ 24 g 
assuming ISM abundances from I Wilms et alj|2000h . 

If v z » v w 1000 km s _1 and i? w ind = 46 lt-s, then 
the implied wind mass loss rate is M win d < 3.9 x 10 20 g 
s -1 , which is roughly 25 times our measured maximum 
mass accretion rate. In iPaper H we estimated M w ind 
based on spherical symmetry arguments and the ioniza- 
tion par ameter. Our estimate here is consistent with 
IPaper J if the covering factor in the spherical approxi- 
mation is of order unity. This may be more appropriate 
for the heartbeat state than for the other Chandra ob- 
servations, since our average HETGS p-state spectrum 
features a P-Cygni profile with roughly equal emission 
and absorption components. The kinetic luminosity in 
the wind is L W md < 10 36 ergs s _1 , which is far less than 
the radiative luminosity. Still, if the wind is as massive as 
implied by our estimates from equation (fTO)) . we should 
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expect it to exert a significant influence on the accretion 
dynamics in GRS 1915+105. 

iShields et all (|1986f ) argued that very massive thermal 
winds can excite long-period oscillations in the accretion 
disk as long as the mass loss rate in the wind is sufficiently 
high (Mwind > 15M) and increases with the accretion 
rate. As the accretion rate rises, the wind becomes more 
massive and begins to drain the disk. This emptying 
quenches the wind, allowing the disk to refill, producing 
a repetitive cycle. Our M w - m d is potentially high enough 
to produce this long instability, with a period 

P « 3400s a- 7 / 9 M 14 / 9 Tj- 4 8 /3 Mr 7 1/3 , (11) 

where M is the black hole mass in units of Mq, Tics is 
the Compton temperature in units of 10 8 K, and Mn is 
the mass accretion rat e in units of 10 17 r s' 1 . 

Previously, X-rav (IWihns et al I2001D and optical 
(jBrocksopp'. Groot. fc Wilms! I200lh observations pro- 
vided strong evidence that this instability may drive 
soft/hard state transitions in the black hole binary 
LMC X-3 , although no wind w as observed directly 
(|Cui et all 120021 : IPage et all 12001 . In GRS 1915+105, 
thermal winds have been suggested as dynamically im- 
porta nt (|Lee et all l2002t iRau. Greiner. fc McColloughl 
|2003l ). but only recently has the contribution of this 
particu lar massive wind to X-ray variabili t y bee n mea- 
sured (jNeilsen fc Led 12001 ILuketic et all l20Toh . For 
our measured parameters in the heartbeat state of GRS 
1915+105, P ~ 1.5 x 10 6 s, or about 17 days, if Ti C8 = 1. 
To compare this result to the available data on the p os- 
cillation, we locate all observations of heartbeats in the 
entire RXTE PCA archival lightcurve of GRS 1915+105. 
Around each pointing, we search for nearby p states in 
windows of various durations. Statistically, we find that 
90% of the pointed observations in a window of width 
12 + 5 days surrounding a p state observation also show 
heartbeats. We conclude that the period of the Shields 
oscillation P is comparable to the typical occupation 
time for the p state, so that the wind-driven instabil- 
ity may actually be responsible for transitions into/out 
of this state. 

Our results therefore indicate that the accretion disk 
wind in GRS 1915+105 may play an integral role in tran- 
sitions between variability classes, effectively acting as a 
gatekeeper or a valve for the external accretion rate, and 
facilitating or inhibiting state transitions. In this con- 
text, it is notable that all published observations of GRS 
1915+105 in a remotely soft or variable state with enough 
sensitivit y to detect highly i onized abs orptio n lines have 
done so feotani et alj|200Ct iLee et al.ll2002HUeda et all 
120091 iNeilsen fc Led 120091 : lUeda et all 120101 ). We note 
that ILee et aLl (|2002f ) actually detected these lines in a 
relatively soft inst ance of the hard \ state, and that we 
showed in iPaper ]] how the wind is sensitive to the frac- 
tional hard X-ray flux, not the specific state. 

Although the duration of our observation is clearly too 
short to probe long-term changes in the accretion flow, it 
provides us with a useful diagnostic of the slower changes 
in the disk wind. Inspection of a series of time-resolved 
spectra (subsets of our observation 3-10 ks in duration) 
indicates that the average ionization level of the wind 
may decrease over the course of the 30 ks Ch andra ob- 
servation (not unlike the changes detected by ILee et all 
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Fig. 14. — The PCA count rate lightcurve for the heartbeat state, 
labeled to correspond to our physical description of the p cycle in 
Section [6] 

l200llSchulz fc Branddl2002l:IUeda et al.H2009h . The ori- 
gin of this change is unclear, and it may have no im- 
mediate impact on state transitions: PCA pointed ob- 
servations show the source to be in the heartbeat state 
on 2001 May 23 (our observation) and May 30, but not 
May 16 or June 5. The sampling is hardly sufficient to 
determine the time of state transitions, but the ASM 
lightcurves suggest that this p state may have persisted 
for at least 3-5 days before and after the time interval 
presented here. 

6. TRACKING THE CYCLE PHYSICS 

For clarity, we present a brief step-by-step summary 
of our model for the physics of the heartbeat state, be- 
ginning at the minimum of the cycle for ease of narra- 
tion. Each step of this narrative is marked on the X-ray 
lightcurve in Figure 1141 Furthermore, each event is in- 
cluded in Table [3] along with its relevant phase interval, 
location in the text, and observational evidence. The 
precise details apply specifically to this observation, but 
in the future, it may be possible to generalize to other 
instances of this oscillation. 

1. (</> ~ 0.3) A wave of excess material, supplied by 
the high external accretion rate, originates near 
25i? g and propagates radially (inwards and out- 
wards). For R < 100 i? g , the disk is always domi- 
nated by radiation pressure. 

2. {<j) = 0.3 — 0.8) Possibly as a result of local Edding- 
ton effects, the disk responds by increasing its inner 
radius (the slow rise) at ^constant temperature. 

3. ((f) = 0.8 — 1.0) During the soft pulse, the disk lumi- 
nosity increases rapidly. When it reaches its max- 
imum, i?i n drops sharply, possibly at a constant 
accretion rate. This suggests a disk instability at 
the end of the density wave. Incidentally, the max- 
imum disk luminosity in this interval is ~ 90% of 
the global Eddington limit for spherical accretion. 

4. ((j) = 0.0 — 0.1) The instability leads to the ejection 
of the material from the inner disk, possibly the 
hot surface layer, which may produce a flash of 
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TABLE 3 
Physics of the p Cycle 



Oycle 
Event 


Phase 
Interval 


Physics 


Location 
(R/R s ) 


Discussed 
in §x 


Observational 
Evidence 


Minimum 


~0.3 


1: Density wave originates, 
propagates in the inner disk 


25 


|5.1.2| 


Strong cyclic variations in -Ldisk an d 
accretion rate; High ^disk/^Edd 


Slow Rise 


0.3-0.8 


2: Local Eddington evolution 


1 5 


I5.1.U5.1.2I 


Ldi s k, ^?im an< l M grow slowly 
at roughly constant T oba 


Soft Pulse 


0.8-1.0 


3: Disk luminosity rises 


1-5 


14. 1.31 15. 1.21 


max(L^; t) ^ Lt?aa & 




0.9-1.0 


3: Disk becomes unstable, 
possibly due to radiation 
pressure 


1 5 


14.1.31 
I5.1.H5.1.2I 


Unstable: i?; n drops, T a ^ B spikes; 
Radiation Pressure: Ldisk S; ^Edd a 
at constant high M (L oc T 4 / 3 ) 


Hard Pulse 


0.0-0.1 
0.05-0.15 


4: Disk ejects material, 
which collides with corona 
5: Density wave subsides, 
disk relaxes 


1 5 

1-5 


14.1. 31 15. 1.31 

14.1.41 
15. 1.1115. 1.21 


^foldi We plummet; Brcmsstrahlung 

normalization spikes 

L disk and T obs drop, R in grows 

at constant low M (L oc T 4 / 3 ) 


Wind Formation 


0.2-1.1 


6: Intense X-ray heating 
launches, ionizes a wind 
from the outer disk 


10 5 - 10° 


[531 


Fexxv, Fexxvi absorption lines 
grow and then fade; Line widths and 
blueshifts are non-relativistic 


Hard X-ray Tail 


0.15-0.36 


7: Production of a 
short-lived jet 


1 - 10 3 


[331 
I4-I.ll 15.1.31 


Band-limited noise, high rms; 

hard X-ray spectrum; low luminosity 



Note. — Physical processes in the p cycle arc arranged here by phase and by location in the accretion disk, along with observational evidence 
and relevant sections of the text. Quoted phase intervals arc approximate. 
a Here LEdd is the global Eddington limit, i.e. the Eddington limit for spherical accretion. 



brcmsstrahlung as it collides with the hot corona 
(the hard pulse). 

5. ((f) — 0.05 — 0.15) The wave of material subsides; 
the disk radius moves out quickly at the new (low) 
accretion rate. 

6. ((f> — 0.2 — 1.1) The bright X-ray pulses travel out 
along the disk, providing an impulse of X-ray heat- 
ing and producing a hot, massive wind, which is 
subsequently over-ionized. Radiation thus links the 
dynamics of the inner and outer accretion disks 
from R=l- 10 6 R s . 

7. ((f) = 0.15 — 0.36) The hard X-ray spectrum, high 
rms, and band- limited noise near hard X-ray tail 
and the minimum of the cycle signal the possible 
production of a short-lived jet. Conditions in the 
accretion disk and corona are similar to other hard- 
state-like dips in GRS 1915+105. 

7. CONCLUSIONS 

In this paper, we have presented a detailed analysis of 
the spectral and timing variability of GRS 1915+105 in 
the 'heartbeat' (p) state, a 50-second quasi-regular oscil- 
lation that usually consists of single- or doubl e-peaked 
burs ts approaching the Eddington luminosity (|TCS97l 
IBOOI ). The accretion dynamics in this state are partic- 
ularly interesting because we find strong periodic vari- 
ations in both the X-ray continuum and in the accre- 
tion disk wind. Our results constitute the very first 
observational probe of accretion disk wind physics on 
single-second timescales, nearly four orders of magnitude 
shorter than the dynamical time at the location of the 
wind. By performing our analysis with respect to cy- 
cle phase (as opposed to time), we gain special insight 
into the physics that drives the oscillation (the radia- 
tion pressure instability and the local Eddington limit) 
and the origin of the accretion disk wind (transient X- 
ray heating). We argue that our observation constitutes 



a snapshot of an evolutionary process around this black 
hole, in which the wind and its variability play an inte- 
gral role in the changing state of the accretion disk on 
timescales from seconds to weeks or months. 

In the context of previous results, we have been able 
to answer two pressing questions about X-ray variability 
in GRS 1915+105: 

1. What is the origin of the hard pulse in the 
lightcurve of the p state? Theoretical models pre- 
dict single-peaked lightcurves (but see JC05), and 
yet single-peaked cycles account for only 40% of 
the observed cycles (Neilsen et al. 2011, in prepa- 
ration) . 

2. How does the wind know about the variability 
in the inner accretion disk? This questio n is ex- 
tremely important given our argument in iPaper J 
that the wind functions as the mechanism for jet 
suppression. 

Our models of the X-ray continuum suggest that the 
hard pulse in the lightcurve is produced by the radia- 
tion pressure-driven ejection of material from the inner 
disk. As the ejected plasma collides with the corona, it 
either produces a burst of bremsstrahlung or becomes so 
Compton-thick that it scatters nearly all the light from 
the inner disk. Simulations (NRM00, JC05) have pre- 
dicted such ejections, but we believe our result consti- 
tutes the first spectral detection of the ejected plasma. 
This type of analysis may be repeated easily for many 
other variability classes, allowing future detailed charac- 
terizations of accretion/ejection physics by X-ray state. 

Of great interest is our ability to diagnose and under- 
stand these X-ray states in the context of inflow and 
outflow physics. It is particularly important to establish 
a clear link between the variability of the inner accretion 
disk and the wind, without which a connection between 
the wind and the jet would seem improbable. In our 
analysis, we find that the fast variability in the X-ray 
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luminosity of the inner disk can actually produce sig- 
nificant structural changes in the outer disk, producing 
a disk wind on the same timescales. Our calculations 
based on equation (fTTj) provide an estimate of (twice) 
the typical timescale for jet suppression if the wind ac- 
tively quenches the jet from the outer disk, but it is clear 
that both fast and slow variations in the accretion rate 
cannot be neglected in this process. Future observation 
of stat es with significant jet and wind activity, l ike the 
state (jMirabel et all [1998 iNeilsen fc Lei [2001 Neilsen 
et al. 2011, in preparation), will reveal further details 
about the precise nature of the wind-jet interaction. 

Finally, we note that instruments like the High Timing 
Resolution Spectrometer planned for the International 
X-ray Observatory (IXO) will revolutio nize fast variabil- 
ity st udies of both emission lines (e.g. iMiller fc Homanl 
2005) and accretion disk winds in X-ray binaries. With 
high signal-to-noise, good spectral resolution, and excel- 
lent time resolution, it may be possible to perform QPO- 
phase-resolved reverberation mapping with disk emis- 
sion lines and wind absorption lines, not only in GRS 
1915+105, but in a large subset of the known Galactic 
black holes. Furthermore, recent XMM-Newton studies 
of Seyfert galaxies have detected variable absorbers dur- 



ing lo ng or quasi-perio dic dips (NGC 1365. iRisaliti et all 
2009; RE J1034+396. IMaitra fc Millerll2010D . The vari- 
ations, which have been attributed to the orbital motion 
of the warm absorber, suggest the exciting possibility of 
a similar database of variability for supermassive black 
hole systems. Future high-resolution spectral variability 
studies with Chandra and missions like IXO will thus al- 
low us to probe even deeper into fundamental accretion 
physics on all mass scales. 
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APPENDIX 
MEASURING PEAK TIMES 
Motivation and Method 

Our primary concern in measuring the peak times for the p cycle is consistency between RXTE and Chandra. For 
RXTE analysis alone, the single-peaked spectral hardness ratios may provide a simpler method for measuring peak 
times. But since Chandra's limited bandpass does not allow optimal hardness ratios, we measure peak times from the 
X-ray count rate. As discussed in Section 13. 1[ we use a representative cycle as a template for cross-correlation with 
the entire lightcurve; the peaks in the resulting cross-correlation correspond to the peaks in the lightcurve. We smooth 
the cross-correlation with a Gaussian of FWHM 3 seconds (much less than the typical width of a correlation peak) 
and measure the time of maxima via parabolic interpolation. The resulting maxima are preliminary peak times for 
the lightcurves. We use them to create the average phase-folded lightcurve, which we then use as a new template for 
cross-correlation. This process mitigates the possible dependence of our results on our initial choice of template. 

PCA versus Chandra 

Given the high effective area of the PCA and the high flux of GRS 1915+105, we are confident in our peak times 
from RXTE. However, because S/N in our grating spectra does not allow use to use only the time intervals where both 
instruments are on-source, we need to be equally confident in our peak times from Chandra. We actually find good 
quantitative agreement between the Chandra and PCA peak times during the overlap interval: the mean difference in 
peak times is AT = —0.98 ± 0.25 s (where the error is the la sample standard deviation for the 273 peaks detected 
by both instruments). Qualitatively, this means that the Chandra peak times precede the RXTE peak times and that 
the noise introduced into the phase ephemeris by using the Chandra peak times, corrected by AT, is much less than 
the cycle period (0.25 s <C 50 s). The fact that AT is non-zero reflects the different instrumental sensitivities of the 
HETGS and the PCA, since the 2-6 keV lightcurves (normalized to mean count rate) have an RMS difference of 0.1%. 
The fact that AT ~ 1 s is a numerical coincidence, since the 2-4.5 keV PCA lightcurve peaks 1.47 ± 0.19 s before the 
full PCA count rate. Physically, the fact that AT is negative indicates that during this cycle, the hard flux peaks after 
the soft flux. 
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